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This paper is the second part of an experimental study on the characteristics of the fluidized state. Part I [1] 
dealt with the critical mass velocity at which fluidization starts. Part II gives an extensive investigation of the 
rate of heat transfer. 


Summary—In a small scale apparatus the heat transfer coefficient to the cylindrical wall of a fluidizing 
bed (“dense phase "’) was measured under strongly varying conditions of gas and solids. The solids used were 
closely sized samples as well as mixtures of different sizes of carborundum, ironoxide (Fe,O,), coke, lead, 
Devarda’s alloy, and fly ash. The gases used were air, mixtures of nitrogen and hydrogen, carbon dioxide, 
argon, town’s gas and methane. 

For a correlation of the measuring results use is made of a system of dimensionless numbers that may 
be considered as a complete representation of all variables which may be expected to influence the value 
of the transfer coefficient. The gas rate is characterized by the reduced mass velocity G/G, where G, is the 
mass velocity at which fluidization starts. A correlation based on this variable proves to be successful for 
Re <5 and can be written as: 
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With the aid of the generalized shape factor B introduced in Part I of this study, the above mentioned 
relationship can be transformed into a correlation restricted to the region 2 < G/G, < 20. This can be 
written as 

Nu 
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Résumé—Dans un appareil de laboratoire, il a été mvsuré le coefficient de transmission de la chaleur 
vers la paroi cylindrique d'un lit fluidifié (phase dense) a des conditions fort variantes quant au gaz et aux 
substances solides. Les substances solides employées étaient tant des échantillons de grosseur exacte que 
des mélanges de carborundum, d’oxyde de fer (Fe,O,), de coke, de plomb, d’alliage Devarda et de cendre 
volante étant différants en grosseur. Les gaz employés étaient lair, les mélanges azote-hydrogéne, le 
bioxyde de carbone, l’argon, le gaz de ville et le méthane, 

Pour une corrélation des résultats de mesure, il a été utilisé un systéme de nombres sans dimensions 
qu'on peut considérer comme une représentation compléte de tous les variables qui influencent la valeur 
du coefficient de transmission de la chaleur. La vitesse du gaz se caractérise par la vitesse de masse réduite 
G/@, dans laquelle G, est la vitesse de masse & laquelle la fluidification commence. Il parait qu'une corré- 
lation basée sur ce variable est exacte pour Re < 5 et s’écrit comme suit: 

Nu 1 918 / oc 0-36 
dans laquelle 
y 0,028 pour 2 < GG, <0. 


Une plus vaste région (1< G/G,< 30) peut étre dominée par 
y (G/G,) = 0,035 (G/G,—1)"*, 


A Vaide du facteur de forme généralisé B, introduit dans la partie I de la présente étude, la relation 
susnommée peut étre modifiée dans une corrélation limitée & la région 2 <G/G, < 20. On peut écrite: 


dans laquelle, pour 

Re <5, (B Re) = 0,58 (B Re)**. 
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INTRODUCTION 
One of the great advantages of a fluidized bed for 
carrying out heterogeneous reactions between gases 
and solids and gas reactions in the presence of a solid 
catalyst is the property that even at a locally high 
rate of heat production the differences in temperature 
throughout the bed are very small. 

This uniform temperature distribution may be 
traced to the rapid exchange of solid particles between 
the various places in the bed and this, combined with 
the great heat capacity of the fluidized solid mass, 
is responsible for a very “‘fast’’ exchange of heat. 

Moreover the marked turbulence of the solid par- 
ticles causes the coefficient of heat transfer to a cool- 
ing or heating surface present in the fluidized bed to 
be several times as large as in the case of a fixed bed 
or an empty tube, all other conditions being the same. 

The present article describes an extensive investiga- 
tion of the value of this heat transfer coefficient (as it oc- 
curs inapparatus on laboratory scale) under widely vary- 
ing conditions of granular material and fluidizing gas. 


Previous Work 

Experiments in this field have been published inter alia 
by Loowinvk [6], Leva, WerntRaus and GruMMER [4] 
and Mickiey and Tricine Like the present 
article the two former publications deal with heat 
transfer under the condition of “dense phase” 
fluidization i.e. at such a rate of gas flow that the 
expansion of the bed is only relatively small. MickLey 
and TrILitic, however, chiefly study the region 
where the density of the bed depends upon the gas 
rate and is up to 20 times as small as the density 
of the fixed bed (dilute phase fluidization). 

LoGwinvk, who measures heat transfer to a water- 
cooled tube positioned axially in the bed, correlates 
the results of his experiments made with fluidized 
aluminium oxide, carborundum, silica gel and “* SiO,” 
of four different particle sizes and air, carbon dioxide 
and helium as fluidizing gas, in the equation : 
S800G° 

a 

Leva etal. study the heat transfer to the cylin- 
drical wall of the bed and give for their measurements 


h = 


* The above mentioned authors give enough details to 
compare their results with our final correlation. Other papers 
on heat transfer in fluidized beds are by L.J. Jotiey, Fuel 1949 
28 114-115; O. Levensprer.and J.8.Waxrox, Chem. Eng. 
News 1949 27 1999; O. P. Acarwat and J. A. Srorrow, Chem. 
and Ind. 1951 278, 321; J. R. Campset. and F. Rumrorp, 
J.S.C.1., 1950 69 373. 

** See nomenclature at the end of the paper. 
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on various kinds of sand and an iron catalyst, using 
air, helium and carbon dioxide as fluidizing medium, 
the following greatly different relationship: 


h = 0-95 6, G15, (2) 


The most striking feature is the large difference 
in the exponent of the mass velocity of the two 
correlations. In addition, the formula given by Leva 
etal. is improbably simple in view of the great number 
of physical constants of solid and gas which may be 
expected to influence the value of the heat transfer 
coefficient. MickLey and TRILLING used glass spheres 
of varying diameters with air as fluidizing medium. 

Thus they could only study the influence of the 
gas velocity and particle size. They correlate their 
results in two formulas, one for heat transfer to the 
cylindrical wall: 

h 0-072 | (3) 
and the other for heat transfer to a tube positioned 
axially in the bed: 

0-028 


(4) 


A comparison of the exponent of mass velocity 
with (1) and (2) cannot readily be made, since the 
formulas of MickLey and TRILLING include the factor 
o, for the density of the fluidized bed, which quantity 
greatly decreases with increasing gas rate in the 
region studied by them. The real exponent of @ in 
eq. (3) therefore is in any case smaller than 0-26 and 
in eq.(4) this value is even negative. 

It seems, therefore, in no way superfluous to 
extend experimental facts while varying the pro- 
perties of both solids and gases as widely as possible. 

The present article deals with heat transfer 
measurements in a fluidized bed of carborundum, 
coke, iron oxide (Fes0,), Devarda’s alloy and lead t+ 
with air, hydrogen, various mixtures of hydrogen 
and nitrogen, argon, carbon dioxide, town gas and 
raw methane as fluidizing gases. 

Of the two first-mentioned solids a great many 
homogeneous sieve fractions varying in size from 
75 to 800u and a few mixtures of various particle 
sizes were investigated. 

All our experiments were confined to dense phase 
fluidization, the variation in the density of the bed 
being sufficiently small as to consider this quantity 
constant for a given powder. 


+ Experiments with sand failed to give reproducible 
results, which was probably caused by electrostatic effects. 
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EXPERIMENTAL 
Apparatus” 
The very simple measuring apparatus (Fig. 1) con- 
sists of a brass cooler measuring 85cm in inner 
diameter and 10 cm in length and fixed between two 
glass cylinders of the same inner diameter and respec- 
tively 10 and 40cm long. The lower side is closed 
with a removable conical brass inlet in which a filter 
of sintered glass has been provided. 

The apparatus is filled, to just above the cooler, 
with powder to be examined and the fluidizing gas is 
supplied through an orifice, by which the flow rate 
is measured, 

An extended brass section at the top side of the 
upper cylinder and having a diameter of 17 cm, 
served as dust collector. 

Just above the filter in the bed an electrical 
heating coil is provided which develops an amount of 
energy q, Which was determined by current and 
voltage measurement. The bed will now adopt such 
a temperature that this energy can be eliminated 
through the cooler. At a cooling water temperature 
of 13 to 15°C it was always possible to regulate 
the heat production in the coil so that the bed 
adopted appr. room temperature, so that heat losses 
to the surroundings were practically negligible. As 
an additional precaution the remainder of the appa- 
ratus was insulated by means of a coating of wadding 
approx. lem thick. 

Another important advantage of the fact that the 
bed was kept at room temperature during the measure- 
ment is that the amount of heat g, added by the 
passing gas, whose inlet temperature ¢, was measured 
in the conical inlet, is only a fraction of g, and can 
be calculated as a small correction from: 


= we, (t, — (5) 


Thus no special care need be taken to measure this 
inlet temperature. 

The bed temperature ¢, was measured with a 
copper constantan thermocouple (wire 0-5 mm thick) 
the junction of which was situated in the middle of 
the bed half way up the cooler. The temperature was 
determined by direct reading from a galvanometer, 
the other junction and a calibrated thermometer 
being placed in a Dewar vessel at room temperature. 

When previously examining the temperature distri- 
bution in the bed by means of a movable thermocouple 
it was found that in the bed zone surrounded by the 
cooler the temperature differed nowhere more than 


bd 


0-1-0-2° C from that of the middle of the bed. At 
a temperature difference of about 10° C between bed 
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Fig. 1. Apparatus for measuring the coefficient of heat 
transfer. 


2 4 6 12 
Fig. 2. Specific heat measurement carborundum 75-110. 


and cooling surface a measurement of the bed tem- 
perature at one place in the middle is therefore quite 
sufficient. 
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In order to ensure a uniform flow of water inside 
the cooler, it was provided with four parallel feed and 
discharge pipes and filled with glass beads with a 
view to good heat transfer. When varying the water 
velocity from | to 12 1 min, the “overall” heat transfer 
coefficient proved to be constant under given con- 
ditions (see Table 1), so that, as the velocity of 12 1 min 
applied for all experiments, the resistance on the side 
of the cooling water could be neglected and the 
temperature of the cooling surface identified with the 
water temperature measured. 

Tn this way the value of the heat transfer coefficient 
was calculated from: 


h d da 


A(th — te) (6) 


where 4 --0-0264 m? represents the cooling surface 
area. 


Table 1. Control of the influence of water velocity. 


Water temperature “overall” h 


Water velocity 


litres min C sec 


inlet outlet 

239 
238 
237 


240 


15-45 
12-9 
13-45 


12-75 


49 

42 12-9 

6-3 13-4 
110 12-75 


The reliability of the measuring procedure is 
demonstrated by the results obtained with a few 
control measurements, where the heat transfer coef- 
ficient calculated from (6) did not vary with abnormal 
deviations in the bed temperature with respect to 
the temperature of the surroundings nor did such 
variations occur when the heating coil was shifted 
from a position just over the filter to just under the 
cooler (Table 2). 


Table 2. 


Control of the influence of the bed temperature and the position of the heating coil 
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obtained. All gases and the greater part of the solids 
are identical to those mentioned in our previous 
paper “Studies on fluidization I.” 


Solids 
All solids used as well as the required data are given 
in Table 3. 

Since the diameter of the particles has only a 
relatively slight influence on the heat transfer, as will 
be discussed later, we have, for the sake of simplicity, 
not used the effective diameter, introduced in Part I 
of these studies, and have based the diameter on 
screen analysis. 

As diameter d of the homogeneous sieve fractions 
we chose the average of the sieve limits. For mixtures 
of different particle sizes we invariably performed a 
sieve analysis, the diameter being calculated from: 


d Sa di (7) 
where x, is the percentage by weight of the fraction 
with diameter d, (see Part 1). 

The density 9, was determined with a pyknometer 
except for fly ash and coke. Of fly ash, which through 
the microscope appears as beautiful transparent 
spheres (see Figure 1 A), the helium and the mercury 
density were determined separately which two values 
proved to be equal with an accuracy within 1%. 
From this it is obvious that the fly ash particles are 
perfectly gas tight hollow spheres. This means that 
for our purposes fly ash may be considered to consist 
of “non porous’’ material of a very small density 
(600 kg 

Coke powder consists of particles whose surfaces 
show relatively large cavities and holes. Never- 
theless both the helium and the mercury density 
proved to be practically independent of particle 
size (see Part I). In Table 3 
the mercury density is given, 


ty ty G h 


14-7 
144 
14-2 


2128 434 
43 «6215 «6437 
36 S574 445 


2443 31-75 O411 
2335 175 0-418 
20-6 19-4 0-420 


coul the helium density being 10°% 
higher. 

The way in which the density 
of the bed, at maximum porosity 
Was determined has been 


just over filter 
just over filter 
just under cooler 


Data of materials used ( solids-gases ) 
In view of the great number of data to be known of 
the solids and gases used, later on a survey will be 
given of the physical constants to be applied in this 
paper and of the way in which these constants were 


fully described in Part I. 

The specific heat c, of the solid was determined for 
every species of material by means of the heat transfer 
apparatus described above. To this end the powder 
was fluidized with air and a normal measurement of 
the heat transfer coefficient was performed in the 
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Table 3 

Solid Sieve limits (1) d 0. G, (air) 
carborundum 75-00 3180 O-O134 
90-100 3180 0-O0178 

110-125 3180 0-026 

150-175 3180 0-054 

175-210 3180 0-069 

210-240 3180 0-084 
75 110 3180 O-O175 
iron oxide 75-110 5180 0-0288 
coke 75-110 1800 00-0146 
125-175 1800 O-0295 

150-210 1800 0-041 

210-240 23 1800 0-067 

240-300 1800 0-100 

300-420 1800 0-161 

420-500 1800 0-225 

S00 1800 0-350 

Dey. alloy 60-110 : 4320 0-02 

lead 50-210 11120 5360 { 0-097 
952 tly ash + graphite 90-125 0-046 


50-210 

75-210 
75 We 
210-240-50% 
iron oxide 50-300 

coke 50-210 


0-O193 
0-025 
O-O104 


carborundum 


0-033 
O-O195 


way already described. Then the heating coil was 
suddenly switched off while the rate of gas flow was 
From this moment the bed tem- 
perature will fall according to an exponential curve 


kept constant. 


and adjust itself at a temperature ¢, just above that of 
the cooling water. This temperature follows from the 
equation 


h A(t, —t,). (8) 


we,(t, 


When the heat capacity of the bed amounts to 
W, and that of the glass wall with glass filter to W, 
it can be easily derived that the fall of temperature 
of the bed as a function of time 7 takes place 
according to: 


t,—t 


hA+ Wey al 
| Wet We 


Cexp: (9) 


By plotting ¢, 
linear curve is obtained, from the slope of which 


t, logarithmically versus time a 


W, can be caleulated. W,, was determined by weighing 
the glass cylinder (under the cooler) together with 
the glass filter, a correction being made for that 
section of the other cylinder (above the cooler) that 
was under the level of the bed. The value of W, 
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was maximum 15% of W,. In Fig. 2 a measurement 
of the specific heat has been exemplified. 

As a guidance Table 3 gives such values of the 
heat conductivity k, of the various solids as are avail- 
able. From our experimental results to be discussed 
later it will be seen that the influence of this quantity 
on the rate of heat transfer may be neglected. 

Finally the critical mass velocity G, of air and the 
generallized shape factor B derived herefrom (see 
Part I) is mentioned. 

Fig. LA gives the microphotograpbs of fly ash 
and lead. For photographs of other materials see 
Part I. 


Gases 

In Table 4 all gases used are given together with the 
various constants. 

The viscosity « of the various gases was deter- 
mined relatively with respect to air (Part I). 

The heat conductivity k was taken from Perry [8}; 
the values for N,-H, mixtures are based on meas- 
urements by Ipps and Hirst [2]. The value of towns 


Vol.I 
No. 2— 1951 
B 
0-66 
0-67 
0-66 
0-62 
0-69 
0-78 
0-66 
0-60 
0-43 
O-47 
0-43 
0-42 
0-42 
0-58 
0-605 
0-265 
u Ww) 3180 1700 641 Is 
110 3180 1663 641 1s 0-62 
v 120 3180 1825 641 Is 1-04 
x 122 5180 2870 632 
A 1s00 720 741 
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gas and crude methane were calculated, according Table 5 
to Lixpsay and Bromiey [5], from the gas analysis Gue CO, CH, 0, CO H, CH, N, 


(Table 5). 
Towns gas Ol OF 645 23-5 40 


determined by analysing the contents of every fresh 
bomb. The values of « and & for these mixtures were 
determined by interpolation for every composition: 
c, and o were caleulated as weight means. 

°C sec 


4 +-+4+— carborundum - air 
+--+ $+ 
+ o 90 Run NP 
— = 44 125 
| 7180-1795 . 
TT 175-210 e 1-12] 


001 00s 01 os 10 
Fig. 3. Coefficient of heat transfer 4 ex mass velocity 
: for various carborundum fractions with air. 
Fig. 1A. Microphotographs of a) flyash 90-125 u, For the various N,-H, mixtures the Pr number 


b) lead 50-2104. For other materials see part I. - : 
deviates fairly considerably from the value of air, 


Table 4 which enabled us to study also the influence of Pr 


Geo (28°) on the rate of heat transfer 


Experimental results 
H,+2%N, O93 O1765 OF108 11-292 0-504 
80% H, + 20% N, O11 6297 3-976 0-436 
5% H, 35% N, 1495 0458) 2.577) 0-427 
45% H, N, 164 0-064 6-671) L764 O45) experimental conditions for earborun- 


In view of the large number of measurements Table 6 
only gives the heat transfer coefficients and the 


25% H, 75% N, 1-71 1-353 0-526 dum-air. 
Air A complete table giving all the experimental results 
Argon 2235 1650 0528) 0-667 be obtained from the authors. 
Town gas 22 1005 60-608) Discussion or ResuLtTs 


In the following development of a correlation of our 


results all measurements for which Re >5 have been 

The volume ratios of the N,-H, mixtures tabulated — omitted since these values proved to deviate systema- 

are rounded off. The real ratio naturally varied more _ tically from the correlation holding for the region of 
or less from the value tabulated and was therefore viscous flow. This point will be discussed later. 
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Fig. 3 gives h as a function of the mass velocity 
G for the various sieve fractions of carborundum with 
air as fluidizing medium. The equivalent course of 
these curves from the starting point of fluidization, 
the value of A ranging from 100 to 400 .J/m? °C sec, 
is demonstrated in Fig. 4, where the same measure- 
ments have been plotted as a function of G/G,. 

Hence it is obvious when trying to develop a correla- 
tion to take into account the gas rate in the ratio GG,, 
the reduced mass velocity. 

This reduced mass velocity may be considered to 
be a measure of the agitation of the bed, which from 
G)G,— 1 (fixed bed) gradually increases. 

For a given value of G/G, the value of h may also 
depend on the following numbers: 


constants of the solid: o,,,, ¢, and d 
constants of the gas: 0, k, ¢, and mu 
and finally the acceleration of gravity g. 


The 9 constants may, in combination with h, be 
reduced to a max. of 6 dimensionless numbers, which 
were chosen as: 

h.d »  00bmge 
= q = 
(10) 


and 
0 Cpe k 


In a complete dimensional analysis of all indepen- 
dent variables the Re-number Gd would also have 
appeared. However, from dimensional considerations 
it is permitted to substitute this Re-number, via the 
dependent variable G,, by G Gy. 

The numbers mentioned in (10) together with the 
reduced mass velocity may be considered as a com- 
plete set of dimensionless combinations of all indepen- 
dent variables which may reasonably be expected to 
influence the value of the heat transfer coefficient. 

From this set of numbers k,k is directly omitted 
since the experimental results do not indicate any 
systematic influence of the thermal conductivity of 
the solid, notwithstanding the fact that its value 
differs greatly for the various solids, This omission 
proves later to be justified when no systematic deva- 
tions, of equal trend with &,, occur from our ultimate 
correlation. 

So the experimental results will have to be cor- 
related in the functional relationship: 


Nu ~ function (a: Pr, OCp ) (11) 


bm’ Com 
* See note page 60. 
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for which we chose the equation: 


= function (G/G,). (12) 


4 
sec 


carborundum-aw 

7S - 90 Run 1-16 
90 - 100 «6 
no- 125 «6€%-6 
150-175 
175- 210 ~ 1-12 
210- 240 1-11 


— 
10 
1 5 10 50 100 


Fig. 4. The measuring results of Fig. 3 plotted vs the 
reduced mass velocity G/@,. 


1000 


4 


carborundum- air 


75 - 90 Run N2 a 1-16 
90- 100 «a 
10-125 
150-175 « 61-10 
175 - 210 
vw 210- 240 


Fig. 5. Dependency on the diameter for various values 
of G/G, derived by interpolation from Fig. 4. 


Particle diameter 


For the measurements on the various carborundum 
fractions with air as the fluidizing gas all numbers on 
the left-hand side of eq.(12) are constant with the 
exception of Nu and gy’, which vary with the diameter. 
In Fig. 5 where these measurements are represented, 
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h has been plotted as a function of d for three different 
values of GG, (3, 5 and 10). The points which were 


104 t 
3 
$10 © 75- 904 Run NP a 
90-10 . .b1-13 
& 110-125 
9150-1795 .@1-0 
175-210 » # @t-2 
103 L tit 
1 5 10 50 100 
Fig. 6. Diameter correlation of the carborundum 
measurements. 


interpolated from Fig. 4 are very satisfactorily 
correlated by parallel lines with a slope of 0-35 


1000 o 
+ +p o 
°C sec an 
h 
s00}__| an = 
A carborundum - gases 
© 985% Run 16-20 
B00%oHg+20% No . 921-3 
46 No . . 41 
. . 942-48 
50 aif «91-15 
~ coal dioxide 
x argon 949-33 
@ town gas 9 61 - Ga} 
@ methane (raw) » «~ 965-68 
—= 
10 
1 5 10 50 100 
Fig. 7. h vs G/G, for different gases with carborundum 
75-110u. 


so that for the exponent of g’ in eq. (12) is found 
y=— 0-45. 

In Fig.6 the correlation of the diameter rela- 
tionship is demonstrated by plotting h/d’™ (instead 
of versus G/G,. 
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Prandtl exponent 
Fig.7 gives, again as a function of G/G,, the transfer 
coefficients measured for the 75-110 u carborundum 
fraction with various gases. 

For determining the exponents x and z in eq. (12) 
we took advantage of the fact that the product oc, 
is practically equal for all diatomic gases so that for 
the gases air, hydrogen and the hydrogen-nitrogen 
mixtures in combination with carborundum the 
number 0¢,/0,,¢, is constant*. 


OS 
y: 
e 985 %oH2+ | 


800 %e Ha + 200 Yo 
619 "Yo Hp + 381 Yo No 
45.6 He + 544 Yo Np 
235 + 765 % No 


— - - 


4 rt OO} 
x“ 

aos 01 05 10 5 10 


Fig. 8. Plot for determining the Pr-exponent x derived 
by interpolation from Fig. 7. 


For the measurements of Fig. 7 with the above- 
mentioned diatomic gases Fig. 8 gives the number 
Nu/Pr'¢g”® as a function of o for a given value 
of GG, and for x= respectively 1, 0-5 and 0. From 
this Fig. the points of which have again been 
interpolated from Fig.7, it appears that only for 
a= 0-5 a linear relationship occurs. The variation in 
the value of Pr is too small to permit a more exact 
definition of this exponent. The value x =0-4 is 
rather satisfactory too. 


Com “s 

For an exact determination of the exponent z of 
Cpm@ for the 75-110 u fraction of coke, carborundum 
and iron oxide in combination with the diatomic 
gases the number Nu/Pr’*g’’* has been plotted in 
Fig. 9 versus 9 for G/G,=5 and 10. The points are 
again very satisfactorily correlated by a system of 


Exponent of the ratios See and 


* This is the reason why in (10) the number 0¢,/0» y,¢, Was 
chosen instead of the quotient ¢,/¢, which is perhaps more 
obvious. 
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parallel straight lines; from the slope of these lines 
it follows that z=0-18. 

The fact that for x=0-5, y= — 0-45 and z= 0-18 
the results for the gases with equal values of oc, are 
well correlated is proved by Fig. 10 where for all 
measuring results with these gases 

Nu ( com 
has been plotted as a function of G/G,. 

Finally we still have to determine the exponent ¢. 

To this end Fig. 11 represents 


© carborund 
4 coke 
iron oxide 


75-110 
75-110 
75-110 


PE} 


Fig. 9. Plot for determining the exponent z of 0» »,/0 
from the results witn diatomic gases (gc, — constant). 


for the various solid-gas system as plotted versus 


for G/G, = 5 or 10 respectively. Because of the 
extreme values of the physical constants of Devarda’s 
alloy and more particularly of lead, the results of 
these solids, which were only available as mixtures 
of different sizes, are also included here, It will later 
be shown more fully that the results on mixtures 
correlate very well with those on samples of closely 
cut sizes. These lines have a slope t= 0-36, so that 
the ultimate correlation reads: 
_ 
All data points for which Re <5 have been plotted 
in Fig. 12 in this way. In view of the extensive 
experimental data the spread proves to be satis- 
factorily small *. 
* The values obtained with Ar, as plotted in Fig. 11 and 


12, are systematically too high, those obtained with CO, 
too low, whereas the other systems of gases and solids give 


a good correlation with the aid of the number _°?. This 


= (13) 
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From the course of the individua] curves in Fig. 3 
and from the correlation curve in Fig. 12 we may 
distinguish between three regions: 

(1) 1<@/G,<2, the transition region from the 
fixed bed (G/G,=1) to complete fluidization. The 
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Fig. 10. Correlation of the gas properties for two atomic 
gases (Qc, = constant). 


x argon 

© 99% Hz 
ar 
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0005 
Fig. 11. Plot for determining the exponent ¢ 
of ms: 


suggests the possibility that also the quotient x of the specific 
heat values at constant pressure and constant volume should 
be included in the correlation. (x4, = 1,67, xcoo, = 1,30, 
diatomic = 1,40.) The spread in the ultimate correlation 
is too large to permit a definite conclusion to be drawn on 
this point. 
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spread in this region is fairly large in consequence 
of the fact that one powder needs a greater value 
of G/Go than another to adopt a complete mobility. 
For GG,=2 this was the case with all solids 
examined. 
(2) 2<@G,<20, the region with constant slope. 


Chemical 
Engineering Science 


so that for G/G,= 'Y* ~75 the velocity of free fall 
is reached. It is clear that an appreciable decrease 
of the bed density must occur much earlier since the 
gas rate between the particles in the bed at the original 
bed density already for GG,=75% porosity ~30 
reaches the velocity of free fall. 

The region 2<@/G,<20 can _ be 


performed on homogenious sieve fractions. 


(3) @/G,> 20, in this region the slope decreases to 
become probably negative for larger values of GG, 
than those applied by us (MickLey and TRILLING [3]. 

The transition at GG, ~20 is undoubtedly related 
with the transition from dense phase to dilute phase 
fluidization. 

A simple calculation shows that this transition 
must actually occur at such a value of G/G,. It must 
be expected that the bed density will decrease greatly 
as soon as the gas rate becomes comparable with the 
velocity of free fall of the particles. 

For smooth spheres the starting point of fluidiza- 
tion is defined by (see Part I, eq. 6): 


Re (14) 


1350" 
For the velocity of free fall Stokes’ law is valid: 


Re 


18 (15) 


co + described by the equation: 
Nu 1 | [ace Nu 1 /@nm\O18/ OC, \0-36 
This equation is represented by 
T | | cI | | the straight line in Fig. 12. 
| | | | S a larger region of the correlation 
| Ai" may be represented by a_ single 
| * equation, reading: 
= 0-035 (G/G,— 
“ In Fig. 12 eq. (17) is represented 
Devarda’s alloy | | by the dotted curve. 
lead | | Th The equality of the exponent 0-45 
] | of and G/G, in (16) suggest that 
- —+_++-+4-- —— + + 
dy | | T should be transferred to the right- 
i © | | l¢ | hand side of the equation resulting 
| in a new variable g’.G@G,, which by 
0 - 100 means of the generalized shape factor 
Fig. 12. G/G,-correlation of all measurements for which Re <5 and B= = ? , introduced in Part I, 


can be directly reduced to BRe. 

In Fig. 13 all measuring results have been cor- 
related on this basis, with the exception of those 
data points for which GG, > 20 resp. < 2. 

As for the various particle sizes the transition for 
G Gy~20 from dense phase to dilute phase fluidization 
oceurs at different values of the Re-number, it is 
obvious that a correlation based on a kind of modified 
Re-number viz. BRe can never include the points 
GG, > 20. 

This also holds for the transition at G/G,=2 so 
that the equation: 


oc, 
Com 


Nu ( 


0-36 
Pr*\ 6 


== 0-58(BRe)°* (18) 
derived from Fig. 13 is only valid for 2<G/G,< 20 
and BRe< 5. 

The nature of the variation of the data points 
with Re>5 is illustrated in Fig.14, where A has 
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been plotted versus Re for those Is) 
measuring series, which include ex- ] 
periments with Re > 5. 

By way of comparison the “nor- 
mal” curve forcarborundum 75-90 
has also been included. The Figure 
clearly shows that the transition 
to the region of decreasing slope, 
which in the case of the normal 
curve occurs for G/G,~20, is more 
and more shifted to the left side 
according as the curve is shifted to 
higher Re-numbers. The transition 
to decreasing slope now occurs for 
Re ~ 5 so that for those curves for 
which G/G,=2 already corresponds 
with Re > 5, the linear upward slope 
for G/G,>2 even does not occur 
at all. From some measurements 
with glass beads (@ ~1 mm) Ah 
shows a marked decline already 
for G/Gy=2. 

A possible explanation of the 
nature of these deviations for Re > 5 
in the G/G, correlation is the follow- 
ing. In the region of complete 
turbulence (Re > 100) Leva [3] finds 
for the pressure drop in a bed of 
spheres: 


iron oxide 


v Devarda’s alloy 


lead 


Run No a,b...s 
carborundum 


a coke 


“Tio 


10 


ts on homogenious sieve fractions for which 2< G/G, < 20. 


+ + He 


5 


10 


O5 


24 1-46 


from which for the starting point 
of fluidization (see Part I) is found: 


Rey = 0-2) (20) 


For the velocity of free fall New- 
ton’s law is valid: 


Re, = 1-7)/g. (21) 


OOS 


+ 


036 
O1 
Fig. 13. B Re-correlation for all 


+ 


From this it follows that in the 
region of turbulence the velocity 
of free fall is reached already for 
G Gg~8-5, whereas in the viscous 
region @/Gy~70 was found. The | 
transition from dense phase to dilute | } 
phase fluidization which for Re<5 O ra 
occurs for G/Go~ 20 may therefore be 
expected already at G/@G,~2-5. This is in accordance The points for- which Re>5 have been included 
with Fig. 14 where for increasing particle sizes this in the BRe correlation of Fig. 13. From this it 
transition occurs at ever decreasing G/G,-values. proves that these data points are now actually 
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correlated although the transition is revealed by a 
decreasing slope. 

Summarizing we may say that the @/G, correlation 
of Fig. 12 is confined to the region of viscous flow 
1000, 
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sieve fractions are less suitable for fluidization than 
mixtures of various particle sizes, heat transfer 
measurements were carried out with some mixtures 
of iron oxide, carborundum and coke. The results 
of these measurements have been com- 


pared in Fig. 15 with the correlation 


sec 


curve of Fig. 12: the correspondence 


=. 


proves to be quite satisfactory. The 


influence of the diameter (d°™) is too 
small for a reliable conclusion to be 


given on the question whether, as in 
the case of the correlation of the 


starting point of fluidization (see Part I) 


also here the definition according to 


eq. (7) of the average diameter of a 


mixture is the most correct. 


With a view to the largest possible 


variation in the properties of the 
solid it seemed worth while to make 


os 


Fig. 14. Character of the deviations for Re > 5. 


(Re <5), whereas BRe correlation of Fig. 13 is 
subject to the restriction 2<G/G,< 20, 

It would not seem excluded, however, that the BRe 
correlation may also be extended for GG ,> 20 if the 


qT 


50-210 Run 1-8 
75-20 i-4 
50°/o - 75-90 


1-4 
x 1-6 

t 1-123 
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50 

Fig. 15. Check on the G/G, correlation by means of 

some special measurements. 


actual density of the fluidizing bed is introduced as 
a variable instead of the density of the fixed bed. A 
further extension of the experimental data is therefore 
required. 


ConTROL OF THE CORRELATION Founp 
The correlation of Fig. 12 is based on more or less 
homogeneous sieve fractions. Since as a rule narrow 


a few measurements with the light 
fly ash. It proved, however, that 
during fluidization of this fly ash, 
serious electrostatic effects oceurred, so that the 
reproducibility of the values measured for the heat 
transfer coefficient was very small*. By mixing fly 
ash with approx. 10% graphite (conductor!) these 
charging phenomena can be avoided. 

A new trouble now experienced was the increasing 
segregation of the fly ash and the graphite at a 
decreasing gas rate. This phenomenon is probably 
responsible for the slight deviations of these results 
from our correlation (Fig. 15). 

Since k, w and ¢, are substantially independent 
of the gas pressure it was possible to control the 
influence of the gas density 0, as it occurs in the 
correlation, by pressure variation and independent of 
the other properties of the gas. The apparatus was 
not suitable for pressures higher than one atmosphere 
so that our measurements were restricted to lower 
pressures. The resulting data points are also included 
in Fig. 15 and satisfy the correlation reasonably. 

Finally Fig. 16 gives the results of other inves- 
tigators as compared with the correlation found bv 
ourselves. Only in Levas paper the values of the 
starting point of fluidization under the experimental 
conditions of the heat transfer measurements are 
mentioned. In Logwinuks experiments. G, was 
determined at room temperature, so that for every 


* See foot note on page 52. 
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measurement this value had to be recalculated for 
the temperature at which the heat transfer was 
determined viz. 100-300° C. Mickiey and 
do not mention a G, value at all. Since these inves- 
tigators only used spheres for their experiments it 
was possible to calculate G, from (14). 

There is no correspondence between the correlation 
found by us and the results obtained by Leva et al. 
as could be expected already in view of the correlat- 
ing eq. (2) used by them. This great deviation may 
exclusively be due to the different experimental pro- 
cedures, 

A better correspondence exists with LogwinuKs 
results though even here the differences are still 
considerable. This may be due to the fact that Loe- 
WINUK measures heat transfer to a tube placed 
axially in the centre of the bed, whereas our measure- 
ments refer to heat transfer to the wall. MickLey 
and Trittixe found a marked difference for these 
two cases. In addition the accuracy of LogwinuKs 
measurements is substantially limited as a result of 
the fact that the correction for the resistance to heat 
flow at the side of the cooling water is considerable. 

Of Mickieys and TRILLINGs 
transfer to the wall only a very limited number could 
be compared with our correlation since the greater 


measurements of 


part of their experiments were made at GG, > 100 


or Re>5. Although their results are approx. 25° 
lower than ours their data points are located more or 
less parallel to our correlation curve. 

When this paper was ready for printing we came 
across an article by Barre, KLAassen and GIsHLer [9], 
who measured heat transfer of a fluidized bed in dense 
phase to a copper tube placed along the axis of the 
bed. The solids used were iron, sand, silicon, Scotehlite 
beads, a cracking catalyst and alumina, the fluidizing 
gas was air. The relation used by the authors to 
correlate their results is rather artificial and far from 
dimensionless. Their results are included in Figure 16. 
Taking account of the uncertainty in the specific 
heat of the solids, which quantity was not mentioned 
by them and of the fact that their particle diameters 
were determined microscopically, which method 
usually leads to higher values than a screening 
the with correlation — is 
reasonably good. 

Finally it was interesting to check the validity of 
our correlation formulas for the calculation of heat 
transfer coefficients in fluidized beds on a technical 


scale. The authors had the opportunity to carry out 


analysis, agreement our 
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a single measurement in a fluidizer on a semi-technical 
scale. The solid used was calcium sulfite powder 
(50-500 4), the fluidizing gas was nitrogen and the 
temperature of the bed 320° C. 

The height of bed was 2-30m and the diameter 
was 30cm. The heat transfer coefficient to a water 
cooled tube of 1m length and 3-5 em diameter was 
determined by measuring the in- and outlet tempera- 
tures of the cooling water. The result for G@G,—25 
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Fig. 16. Comparison between the correlation found 


and the results of others. 


was 315 Jim? °C see which was only 10% lower than 
the value caleulated from our correlation. Though the 
agreement of a single measurement is quite insuffi- 
cient to establish the accuracy of our correlation for 
the design of large scale equipment, it will certainly 
satisfy the need of a rather accurate estimation. 


NOTATION 


A = surface area of cooler, m*( = 10-76 ft*) 
generalized shape factor (dimensionless) 

(} = mass velocity based on empty tube, kg/m* see 
( = 738 lb ft? hr) 

Gg = critical mass velocity 


G/Gy = reduced mass velocity (dimensionless) 
Nu = . = Nusselt number (dimensionless) 
Pr = Prandtl number (dimensionless) 

Re, Reg = “ resp. “ Reynolds number (dimensionless) 
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Cy ¢, ~ specific heat of gas (at constant pressure) of solid 
J/kg °C (0-00024 BTU Ib F) 
d- particle diameter based on screen analysis m 
( 3-28 ft) 
q~ acceleration of gravity 9-8 m/sec? (4-15 » 
10 ft/hr®) 
h heat transfer coefficient °C sec 
0-176 BTU/hr ft® F) 
k. k, — heat conductivity of gas or solid J/m*°C see 
0-578 BTU ft ° Fh) 
Ge % «quantity of heat per unity of time delivered by 
electric coil resp. by flowing gas J/sec 
3-415 BTU/h) 
tir ty t, temperature of gas, resp. of bed, or of wall © 
w mass flow rate kg/sec ( — 7920 Ib/h) 


45 porosity of the bed (dimensionless) 
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pressure drop through bed N/m? ( = 8-7 
10° |b force, ft*) 

pw — dyn. viscosity of gas kg/m sec ( 
density of gas or of solid kg/m* 


2420 Ib/hr ft) 


( = 0-0624 Ib /ft*) 
Chm ~ density of bed at maximal porosity 


0 density of fluidizing bed 
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Estimation of the free enthalpy (Gibbs free energy) of formation of organie compounds 
from group contributions 
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Summary—A method is proposed whereby the free enthalpy of formation of organic compounds in the 
ideal gaseous state is resolved into contributions attributable to atomic groups. 

Free enthalpies of formation of complex molecules are readily estimated by summation of the contrib- 
utions of their component groups. The average deviation of the free enthalpy of formation of organic 


compounds caleulated by this method is 


I. Ixrropverion 
Now and again every chemist and chemical engineer 
will face the problem whether the preparation of a 
particular product according to a new reaction- 
principle is possible. 

Chemical thermodynamics, enabling chemical equi- 
libria to be caleulated, can give the solution to this 
problem, provided sufficient data be known about 
reactants and reaction products. 

In fact, the thermodynamic potentials constitute 
the driving forces causing every natural process to 
proceed in the direction of its eventual state of equi- 
librium. 

Whether or not a reaction will proceed depends 
on kinetic factors. A certain amount of activation 
energy is necessary to keep up practically any reaction. 
However, in many of the cases in which a reaction 
was thermodynamically possible it also proved 
possible to find a catalyst that was sufficiently 
active and selective to realize this reaction. 


06 keal mol. 


The combined knowledge of the thermodynamic 
driving forces and of the corresponding rate factors 
presents the stock from which the chemical engineer 
may take the data enabling him to calculate con- 
verters. Thermodynamics determine the possibility, 
kinetics the rate of the process. 


Unfortunately the application of thermodynamics 
is in many cases handicapped by a lack of sufficient 
data. In these cases it is important to have at hand 
a simple and quick method for calculating these data, 
if only by approximation. This will provide an 
insight into the ranges of temperature and pressure 
in which to look for a catalyst. 


Il. Taermopynamic 
AND EQUILIBRIUM 


Given the reaction of formation of a compound (C 
from its composing elements E, and E, 


nk, 


mE, 


(II-1) 


IL. 
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The equilibrium of reaction (II-1) at a given tem- 
perature, is given by the equilibrium-constant 
(C) 
Thermodynamics gives the relation between the 
equilibrium-constant and the free enthalpy of for- 
mation 


(I1-2) 


AG, =—R.T.InK,. 


From (11-2) we see that the equilibrium in (IT-1) is 
shifted to the right if A,> 1 and to the left if K,<1, 
from which it follows that—bearing in mind (II-3)— 
the equilibrium lies to the right in (II-1) if 4G, is 
negative and to the left if AG, is positive. In other 
words: at a given temperature a compound is stable 
if at that temperature 41G,<0. 

The value of AG, is a measure of the stability of 
a compound with respect to the elements. If it be 
asked to define the position of the equilibrium of the 
reaction 


(11-3) 


C+DeF+G6 
this equilibrium is determined by the constant 


- (F).(G) 


(I1-4) 


(11-5) 


Which again bears relation to the change in free 
enthalpy in consequence of the reaction 


AGreacton = —R.T (11-6) 
in which 
A Greaction =A + A {4 Gy, + A Gio} (IL-7) 


If the free enthalpy of formation of the com- 
pounds participating in the equilibrium is known, it 
is possible to calculate the position of this equilibrium. 

From the preceding it follows that it is of great 
importance to know the numerical value of the free 
enthalpy of formation. 


III. Structure anp THERMODYNAMIC 

PROPERTIES OF ORGANIC COMPOUNDS 
Only a very small part of the overwhelming number 
of known organic compounds have been examined for 
their thermodynamic behaviour; besides it is highly 
improbable that the thermodynamic data of all these 
compounds will ever be measured exactly. Hence it 
is obvious that methods have been and are being 
sought to calculate these data. 

Theoretically it is possible to calculate these data 
according to a statistical-mechanical method. How- 
ever, this method is very laborious and moreover the 
spectroscopic data required to this end are usually 
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lacking. Of late years Rossrx1 and co-workers [39] 
have performed excellent work on the subject in the 
National Bureau of Standards (A.P.I. Research 
project 44). 

For years a method has been sought to calculate 
the numerical values of thermodynamic properties in 
a simple manner and with the help of as few data as 
possible. More or less rough correlations between 
structure and thermodynamic properties have been 
known for a long time, as may be seen e.g. from the 
work of Karascu [30] and Parks and Hurrman [34]. 
Only of late years did it prove possible to state a more 
or less accurate relation between structure and 
thermodynamic properties of an organic compound. 
Some methods will be mentioned here: 


1. The method of Anderson, Beyer and Watson {1\ 


These investigators consider every compound to be built up 
of a base group in which cther molecular groups were sub- 
stituted for the atoms in the base group. They make correc- 
tions for contributions from double bonds, cis- and trans- 
isomerization, conjugation effects and ortho-meta and para- 
substitution. In this way they consider ¢.g. all the paraffins 
as derived from the base group methane by substituting the 
hydrogen atoms by CH, groups. In a similar manner they 
obtain values for the secondary amines by starting from the 
base group HN(CHs,),. 

The disadvantage of this method is that it only gives 
values for the heats of formation and the entropies at 25° C, 
so that the free enthalpies of formation, the heats of formation 
and the entropies at other temperatures have to be calculated 
with the help of a correlation of specific heat given by them. 
Moreover the result of their method depends upon the way 
selected in building up the equation for the specific heat. 


2. The method of Souders, Matthews and Hurd {46} 


Sovupers, Marruews and Hurp consider specific heat, heat 
content, heat of formation and entropy of a hydrocarbon as 
built up of vibration-, rotation- and translation-contributions 
of the composing molecular groups. In this way a thermo- 
dynamic magnitude of 2-methyl-pentane is ¢.g. built up of: 
(a) Vibration contributions from: 
3 —CH, groups, 2 —CH, groups and | —CH group 
(b) Internal rotation contributions from: 
3 groups rotating about the —CH group, | —-CH,—CH,— 
rotation and 1 CH,—CH,— rotation. 
(c) Translation contribution amounting to 4 R. 
(d) In calculating the entropy a symmetry correction amount- 
ing to Ring (see note page 69). 
This method gives very accurate results (see the original 
publication); however, it has the disadvantage of requiring 
a rather extensive knowledge of the molecular structure. 
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3. The method of Bremner 
and Thomas {9} 


Bremxer and Tuomas developed a method — Temperature 
for calculating the free enthalpy of formation x 
of organic products based on the assumption 

that the molecule is built up of characteristic 

groups each of which giving a certain con- 208.16 
tribution to the free enthalpy of formation. pom 
Corrections for conjugation and mutual ase 
influence of various groups are not given ene 
by these investigators. However, they do 

give corrections for the formation of cyclo- =e 
hexane, cyclohexene, the furane ring and 1000 
the tetrahydro-furane ring. The comparison sane 
between the free enthalpy of formation 4 


calculated according to their method and 


H®) cal. mol.» 


i C,H,OH cH,—C=0 CH,COCH, 

cale. lit. (44) | [21] 

51-50 —52.85 | —58-69 —57-53 

64-82 64-68 55-03 56.03 62-64 61.67 
68-64 —68-81 —58.04 58-67 

60.62 61.05 69.22 68.55 

6514 65.28 74.99 74.40 

69-12 69.038 80.19 79-72 

85-44 84-61 

$9.95 88-93 

92.17 


the known enthalpy is very concise, so that 
it is impossible to judge the value of this method from 
their publication. However, a comparison of this method with 
the other methods of calculation mentioned here renders 
the accuracy of the Bremner and THomas correlation 
doubtful. 


4. Rossini’s rule 


A rather accurate correlation for calculating thermodynamic 
magnitudes is applied by Rosstnt and co-workers in compiling 
the A.P.1.-N.B.S. tables for the thermodynamic magnitudes of 
gaseous hydrocarbons (see references given in [39]). These 
investigators consider the thermodynamic magnitude as built 
up of the corresponding magnitudes of compounds with 
analogous structure, ¢.g. 


CH, 
CH,—C—CH—CH, 
CH, cH, (1-1) 
CH,—C—< ‘H, CH, 
CH, 
in which @ is equal to Ine; R. Ina; 


CY; 1H) or 
Rosstxi and co-workers found that this rule gave excellent 
results with hydrocarbons. 

It is of importance to know whether this rule is also 
applicable to non-hydrocarbons. Prrzer [6] states that it is 
possible to calculate the thermodynamic magnitudes of the 
higher normal mercaptans in this way, starting from ethyl 
mercaptan. In this manner the thermodynamic magnitudes 
of three compounds were calculated, by which it was found 
that the magnitudes obtained in this way tally very well 
with the corresponding magnitudes calculated from spectro- 
scopic data *. 

* Calculated were the freé enthalpy functions of isopropy! 
aleohol, acetic aldehyde and acetone: 
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Free enthalpy functions were chosen in order to prevent 
the heat of formation from being introduced into the caleu- 
lations, for as a rule the heat of formation is known only with 
insufficient accuracy. Due to a lack of necessary data it prov- 
ed impossible to carry out more extensive control calculations. 
The calculations made, however, are sufficient proof of the 
validity of Rosstnis rule for the carbon-carbon bond. 


From this it may be concluded that this method, in as far as 
carbon-carbon bonds are concerned, gives rather accurate 
results. Prrzer (6) proposes a calculation of the thermo- 
dynamic data of dimethy! sulphoxide and dimethyl! sulphone 
according to: 


So, SO, 


—® 


4 
(CH,),80 CH,—S—CH, 


(111-2) 


=@® 
(CH,),80, (CH,),SO SO, 80, 


Owing to a lack of data it is impossible to check the validity 
of Rosstxts rule for bonds other than carbon-carbon. 


5. Franklin's method (15) 


FRANKLIN considers the molecule to be built up of molecular 
groups, each of which giving a particular contribution to the 
thermodynamic property. He gives group contributions to 
the heat of formation, heat content, free enthalpy function 
and free enthalpy of formation. In calculating the free ent- 
halpy function and the free enthalpy of formation he applies 
a correction for symmetry (see note page 69). 

For hydrocarbons FrRaNKLINs method gives results within 
2 kcal/mol. The value of the few group contributions he gives 


CH,—C—OH CH,—C—H C,H,OH 
CH, CH, 
H 
Oo 


CH,—C—CH, CH,-C-O 
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for non-hydrocarbons is doubtful. A more extensive discus- 
sion of FRANKLINs method, which is a very simple one, will 
be found hereafter. 


IV. Mopiriep FRANKLIN CORRELATION FOR THE 
FREE ENTHALPY OF FORMATION OF ORGANIC 
COMPOUNDS IN THE IDEAL GasEous STATE 


Starting from a statistical-mechanical calculation by 
Prrzer (35) FRANKLIN proved that a thermodynamic 
function of an arbitrary paraffin hydrocarbon may be 
built up as the sum of contributions from characte- 
ristic groups, of which that hydrocarbon may be 
conceived to be built up, and, if required, with the 
addition of a correction for the symmetry of the 
molecule. As this method gave excellent results for 
the paraffins FRANKLIN assumed that it also held 
good for the other hydrocarbons, which afterwards 
appeared to be true. So for the free enthalpy of 
formation of a hydrocarbon we may put: 


4G, = ine (V-1) 


Composing groups 


in which R.T . Ino represents the correction for the 
svmmetry*. The equation (IV-1) did not give 
accurate results over the entire range; for certain 
forms of molecules constant deviations were found. 
In order to make up for this FRANKLIN introduced 
corrections for these forms of molecules, so that 
consequently (IV-1) changes into: 


AG oF contributions of 14 v corrections 
‘ne composing groups necessary 
+R.T.Ino. 


FRANKLIN denotes contributions to the enthalpy of 
formation at 0; 298-1; 400; 500; 600; 800; 1000; 
1200 and 1500° K, which he calculates from the first 
terms of the homologous series. 

FRANKLINs method gives information only for the 
temperatures at which contributions are given. In the 
present paper this is obviated by considering the group 
contributions as a linear function of the temperature. 


(IV-3) 


AG, =A + 100 . . 


group 


* The symmetry number. 

This is: “The number of indistinguishable positions in 
space the molecule may be made to take up by a simple rigid 
rotation.” 

The symmetry number of a molecule may be found by 
considering a space model of the molecule concerned. A 
great number of symmetry numbers is given in HerzBERGs 
publication [25}. 
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This was based on an argumentation by SCHEFFER [43 }. 
The equation (IV-3) shows a strong similarity to the 
general equation 


AG=AH—T.AS. (IV-4) 


If ([V-3) and (IV-4) are compared it follows that A 
has the dimension of a heat of formation and B/100 has 
the dimension of an entropy of formation. According 
to Uxicn [54] A is by approximation equal to the 
heat of formation at 298-1° K and B/100 approximately 
the entropy of formation at 298-1° K. 

It proved impossible to describe the variation of 
the contribution to the free enthalpy of formation in 
the entire temperature interval (300 =1500° K) by 
one linear temperature function. However, it was 
possible to represent the temperature dependence for 
all the contributions by two straight lines in the 
ranges 300 — 600° K and 600 — 1500° K. 


V. Mernop or CALCULATION 


The values for the free enthalpy of formation are as 
a rule stated for the ideal gaseous state at a fugacity 
of l atm. This also holds for the group contributions 
given in this paper. 

As a basis for the calculation of the free enthalpy 
of formation contributions for hydrocarbon groups 
the free enthalpies of formation as given in the 
publications of the A.P.I. research project 44 {39}, 
[40], [41] and [42] were chosen, with the addition of 
data concerning the free enthalpies of formation of 
cyclopropane [2] and naphthalene [6a]. 

Group contributions for non-hydrocarbon groups 
were calculated from data found in various publi- 
cations. For the non-hydrocarbons only the free 
enthalpy function and the heat content are given in 
most cases, from which the free enthalpy of formation 
is calculated in the following manner: 


AG, = 7.A( + — (V-1) 


in which 


A(H®— Hg) 


elements 


The functions (= and (H®= are as a rule 
calculated from spectroscopic data and are fairly 
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accurate. The magnitude A H) is calculated from the 
heat of combustion. In calculating the group con- 
tributions no values for 4H) or for the heat of 
combustion of the compound were found in some 
cases. Then the heat of formation was calculated by 


Table 1. Alkane-groups 


Group contribution A + 


100 
keal mol 


300 600° K 600 1500 K 


B 


1300 1$00 
Temp °K ——> 


700 900 1100 


Fig. 1. 


means of the correlation given by ANDERSON, Beyer 
and Watson [1], which for this magnitude corresponds 
rather well with the real values. Furthermore for the 
non-hydrocarbons, data were not always available 
up to the 1500° K limit. It was assumed that in these 
cases the temperature dependency up to 1500° K is 
approximated with sufficient accuracy by the formula 
given. In these cases the highest limit at which data 
from the literature were still available is denoted. 


VI. Tue Group Contrrisutions 


1. Hydrocarbon groups 


Contributions to the free enthalpy of furmation were 
calculated for the various groups. In as far as possible 
average values were chosen for the contributions. 
The separate classes of group contributions will be 
discussed hereafter. 


(a) Alkanes 

As starting point for calculating the group contri- 
butions we chose the value for the —CH,— contribution, 
as given in [39]. Starting from this we calculated, 


Table 2. Alkene-groups 


Group contribution A +- ~ > 


100 
mol 


300 600° K 
A B 


+ 11-552 
13-737 
16-467 
17-663 
17-187 
20-217 
25-135 
45-250 
49-377 
51-084 
52-460 


scect 


2 


2 


900 noo 1300 
Temp °K 


300 
Fig. 2. 


by means of the formula (IV-2), contributions for the 
—CH,, —CH and —C—groups. The results of these 
calculations are represented in Table 1. 


(b) Alkenes 

Starting from the group contributions represented in 
Table 1 we calculated group contributions for alkene 
groups. It proved not possible to calculate contri- 
butions of sufficient constancy for the H,C = group. 
However, it was possible to calculate contributions for 
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Group 
A B 
Group 
A | 4 B 
H,C CH, 145 9450-888 
—CH, 10-943 12310-2436 H.C = 1.915 16-255 1-966 
—CH,—- | — 5193 2430 5-830 2-44 Ho 1.965 16-180 2-116 
| 0-705 2.910 0-705 2.910 1915 158152082 
—C— + 1958 3-735 4 4385 3-350 2-205 | 19-584 2-354 
>C=C< 2-573 25-1352. 73 
VOL, 
60 1.027 43-634 l 
1.035 48-170 1.208 1952 
ao} | HC =C=C< 1474 
1-483 52-460 1.483 
60 
HCH 
60 
— = 
20 
q 
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the -+C*~ group from a free enthalpy function of 


naphthaline [6a] and the heat of combustion of 
naphthalene [28]. 


the groups H,c—CH; H,c—c<; Ho—cH; 
>C=C<; H,C=C=CH; H,C=C=C< and 
Ho—c—cH. 


(c) Conjugated alkenes 
In order to avoid correlations for conjugation effects, 
which correlations are neither constant, contributions 
were calculated for characteristic resonance groups 
that will be denoted by a double-headed arrow. 
Contributions were calculated for the H,C--; 
and >C+ > groups. 

Table 3. Conjugated alkene-groups 

Group contribution A 

keal/mol 


300 — 600° K 600 — 1500° K 
Group — 
A B A B 


0-675 
1-035 


4-500 
6-980 


Ho.. 


1-505 | 10-370 


(d) Alkynes 

As starting point for the contributions for the alkyne 
groups we chose the half of the value of the free ent- 
halpy of formation of acetylene after subtraction of 
the correction for the symmetry. With the help of 
the contribution for the HC= group obtained in this 
way we calculated the contribution to the free ent- 
halpy of formation for the —C= group. 


Table 4. Alkyne-groups 


Group contribution A 
keal/mol 
300 ~ 600° K 600 ~ 1500° K 
Group 
A B A B 


HC= 


27-048 
26-938 


0-765 26-700 
0-525 26-555 


0-704 
0-550 


(e) Aromatics 

As was the case with the conjugated alkenes, we also 
calculated contributions for the resonance groups for 
the aromatics. HC3~ was calculated from benzene 
and with the help of this group the —C7— group was 
calculated. For the condensed aromatics we calculated 


Table 5. Aromatic-groups 


Group contribution A 4+- 


T 


600 1500° K 


HC 
Cz 


2-505 
5-010 
6-260 


({) Ring formation 

In order to create the possibility of calculating the 
free enthalpy of formation of the cyclo-alkanes and 
cyclo-alkenes it was necessary to calculate corrections 
for the formation of the various rings. Table 6 gives 
the corrections for the formation of the cyclo-alkane 
rings. The value for the formation of the 4-ring has 
been calculated under the assumption that the contri- 
bution for the formation of the cyclo-alkanes is a 
linear function of the angle of deviation (difference 
between the valence angle in the ring and the valence 
angle in the tetrahedron). Table 7 gives corrections 
for the formation of pentene and hexene rings. 


Table 6. Ring formation I 


600 — 1500° K 


600 =~ 1500° K 


A B 

2.395 3-915 —2.250 

2.045 —8-810 —2.071 
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- 
mo 
300 600° K 
A B A B 
3-047 0-615 0-706 
4-675 1-150 0-988 
6-608 0-514 0.583 
OL. 
L 5-437 0-832 
952 7-407 1-088 
9-152 1.308 
B 
Correction A +- 100 = 
kcal/mol 
300 600° K 
oe 
| 
3 ring 23-458 3-045 22-915 | — 2-966 
4 ring 10-73 2-65 10-60 — 2-50 
. 5 ring 4-275 — 2-350 2-665 —2-182 
6 ring 1-128 1.930 | — 1-504 
Table 7. Ring formation I 
Correction A dl my 
100 
keal/mol 
300 600° K 
Group | — 
A 
Pentene ring ~ 3-657 
Hexene ring ~— 9-102 
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Table 10. Branching in aromatics 


(g) Branching 


In order to obtain sufficiently accurate results, it Conestin A + Bor 

proved to be necessary to calculate corrections taking initia 

into account the branching influence. In this way anu | aapenaea 

corrections were calculated for the influence of branch- Group ; 
A B A B 


ing in paraffin chains (Table 8), for the influence of 


branching of cyclo-paraffins (Table 9) and in aromatics 
(Table 10). Double branching 


Table 8. Branching in paraffin chains 


300 — 600° K 600 — 1500 K 
Groups 1,3,5 position 0 0 0 0 
A B A B 


Side chain with 2 or more 2. Non-hydrocarbon groups H 
C-atoms. 31 Contributions were calculated for the —OH; —C=O; 
3 Adjacent —CH— groups. 2-12 0 212 0 O O O 


H 
—C—; C=C=0; >C=C=0; —C—OH; —C—O—; 


Adjacent —CH and 
—O—; —C=N; —N=-C; —NH,; >NH; —N<; 
« 1-80 0 1-80 0 —NO,; —SH; —S—; >SO; >SO,; sx; 
3 Adjacent —C— groups | —F; Cl; —B, and the —I groups. Furthermore a 
correction was calculated for the formation of 
Table 9. Branching in cyclo-paraffins he and > rings. The free 
Correction A + | b< 
mtn ~ enthalpies of formation were considered as being built 
300 - 600° KK | 600 = 1300°K up of the contribution for hydrocarbon groups with 
Grow p the addition of a contribution for the non-hydrocarbon 
4 group and, if necessary with a contribution as a result 
Branching in 5 ring of symmetry. 
Single branching — 144 0 1-69 0 (a) Oxygen compounds 
@ In contradistinetion to the ideas of Frawx- 
12 positi es | e Lin [15] and Parks and Hurrman [34] only one 
tans 1,2 position |-2:55 | 0 0.945 0.296 contribution for the —-OH is given. Neither are there 
cis 1,3 position 120 (OO _ 0.370 — ogg corrections given for secondary, tertiary and aromatic 
trans 1,3 position 2.35 0 0.800 60.264 —OH groups. This view is upheld by the validity 


of Rosstnts rule (in the correlation proposed here the 
Branching in 6 ring place effect is calculated in the hydrocarbon group). 
For it is possible to calculate the thermodynamic 


Single branching 0.93 0 0-230 —0-192 

Double branching magnitude of isopropyl alcohol according to 
1,1 position 0-835 0.367 1-745 — 0.556 CH, 
trans 1,2 position 241 0-045 0.398 CH, CH, 
cis 1,3 position 1647 0-185 
trans 1,3 position io 8 0.260 0.290 Additional contributions in consequence of place 
cis 1,4 position 0 effects as given by and Parks and Hurr- 
trans 1,4 position -—280 © 0-095 —0-245 are to be attributed to errors in the determination 
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1,3 position O31 031 0 
eee 1,4 position 0-93 0 0-93 0 
Correction A + 100 m4 Triple branching 

keal mol 1,2,3 position 1-91 0 2-10 0 
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of the heat of combustion. The contribu- Table 11. Oxygen containing groups 
tion for the —OH group was calculated un bution A+. Boer 
from ethanol [44], [61]. Methanol shows keal/mol 100 
a deviation, as appears from the follow- 
ing series of free enthalpies of forma- 
tion at 25°C, calculated from data by 
WaLkKeEr [58). 
CH,OH AG} = —38-7 keal./mol. 
C,H,OH AG} = —40-3 keal /mol. 
nC,H,OH AG} = —39-9 keal./mol. 
nC,H,OH AG} = —37-5 keal./mol. 
n€,H,,0H AG? —35-5 keal./mol. 
nC,H,,0H AG} = —33-7 keal./mol. 


H 

—C=—O This contribution was caleu- 
lated from the free enthalpy of forma- 
tion of acetaldehyde (36). 

>C=O Calculated from the free 
enthalpy of formation of acetone given by 
Scuuman and Aston [44], and Gopnev, 
Payvukuina and [21]. 92-62 

Ho—C—0 and >C—C=0 By apply- 14-515 
ing Rossinis rule on the free enthalpy 12°86 
function of H,C=C=O [13] the cor- | 9-462 
responding functions for 


and (CH ,).=C=C=0 were calculated. 


| Highest 
"temperature 
300--600°K 600 ~1500° K of lit. value 


H 
* “cu,—c-cu, °cH,—cH, 


(VI-2) 


—C—OH This contribution was calculated from 
the free enthalpy function of acetic acid, which was 
again calculated from the corresponding function of 
formic acid by means of Rosstnts rule [62]. 


"cH, COOH * “cH, CHO (VE-3) 300 S00 700 900s 100 1500 
O Temp °K ——> 
ig. 3. 


—C—O— Vvepensky, Yvannikov and Nekra- 
sova [56] 1 the position of the equilibrium of =—2— VaLawris [55] determined the positon of 
the equilibrium 
2 C,H,OHyy) C,H,OC,H, (9) +  (VI-5) 
+2Hy- (VI-4) from this calculated the change in free enthalpy 
The contribution from this group was calculated of this reaction. From this magnitude the contri- 
from data by these investigators. bution for the —O— group was calculated. 
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3 °K 
16 1500 
8 1500 
600 
4) 
54 1500 
| 0-77 | 30-15 0.83 1000 
0-91 —28-08 0-91 1500 
| 
2.473 90.569 2.958 1500 
2-86 98-83 1500 
OL. | | = 
| 
P| 2-61 | — 92-62 2-61 800 
0-295 -— 14-515 0-295 900 
0-46 | — 12-86 0-46 900 
0-72 |— 938 0.73 900 
20 
. 
( CH, 2 -60 
|_| 
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Table 12. Nitrogen containing groups O>=~ Bremner and Tuomas [9] give 
= the heat of formation and the freeenthalpy 
healimel of formation of furane. Assuming Uticus 

rule [54] to be valid in this case, the 
contribution from the Os" group was 
calculated with the help of these data. 


300 600° K 


(b) Nitrogen compounds 

and —N_C These groups 
were calculated respectively from the 
free enthalpy functions of methyl cyan- 
ide [14], [35] and aerylonitril [27] and 
from the free enthalpy function of methyl 
isocyanide [14]. 

—NH,; >NH and >N— Contribu- 
tions for these groups were calculated 
from the heat of formation given by 
Bicnowsky and Rossryi and entropies given by 
Aston and co-workers [3], [4], [5], after which the 
free enthalpies of formation obtained in this manner 
were extended to higher temperatures by means of a 
specific heat correlation given by Dopratz [12]*. 

— NO, Calculated from the heat of formation 
given by Bichowsky and Rosstn1 [7] and the entropy 
given by Prrzer and Gwryw [36]. Contributions for 
higher temperatures were calculated with Dospratzs 
specific heat correlation. (The C) equation obtained 
in this way was adapted to some C) values measured 
by Prrzer and Gwryy [36].) 

N= The contribution for this group was calculated 
from the free enthalpy function of pyridine [31]. 


(c) Sulphur compounds 

—SH_and —S— Contributions calculated respec- 
tively from the free enthalpy functions of methyl 
mercaptan [6], [8] and ethyl mercaptan and from the 
free enthalpy function of dimethyl sulphide. 

“8 The contribution to the free enthalpy of 
formation of the resonance group >~S was calculated 
from the free enthalpy function of thiophene [57]. 


* Dospratzs correlation for the specific heat is based on 
the, approximately right, assumption that a given bond 
between two atoms shows two characteristic frequencies of 
vibration, which are independent of the nature of the molecule 
in which that bond occurs. In his publication Dopratz gives 
values for these frequencies for the various bonds and contri- 
butions of these frequencies to specific heat. By the further 
assumption that every degree of internal rotation contributes 
the full value of R/2to the specific heat this investigator arrives 
at a specific heat correlation for gaseous compounds that 
corresponds rather well with the reality in the case of not 
too complex molecules. An extensive explanation of this 
correlation may be found in [12] and [26]. 


HCN 31-179 0-826 30-874 0-775 
—<C=N 30-75 0.72 30-75 0.72 
—N=C 46-32 0-89 46-32 0-89 
NH, 11-606 2-556 -12-972 2-784 
. -NH, 2-82 2-71 - 6-78 3-98 
NH 12-93 3-16 12-93 3-16 
>N— | 19-46 3-82 19-46 3-82 
| 11-32 1-11 12.26 0.96 
—NO, 9.0 3-70 14.19 4-38 
80 
mat VOL. 
—* 
| 
— 
NH 
“20 — -0n 
-40 
300 500 700 1800 
Temp. °K ——> 
Fig. 4. 
80 
40 
7—~CEN 
HCN 
HC= 
30 500 700 «900 
Temp °K ——> 
Fig. 5. 
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>8O and >SO, Contributions for these 
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Table 14. Halogen containing groups 


groups were calculated from the free enthalpy 
functions of dimethyl sulphoxide and di- 
methyl sulphone, which in their turn were 
calculated according to (III-2), following a 
suggestion by Barrow and Prrzer [6}. 


(d) Halogen compounds 

=F Tuompson and Tempe [53] give 
the function of the free enthalpy of 1,1,1 
trifluorethane, from which the contribution 
for this group was calculated. 

=Cl The contribution for the —Cl group 
was calculated as the mean value of the 
increments from the free enthalpy functions 


Group contribution A T 


100 


1500 

1200 

0-397 1500 
15-26 1000 


of methyl chloride, methyl dichloride, methyl 
trichloride, carbon tetrachloride [49], dichlorethane 
[24], vinyl chloride, cis-dichlorethylene, trans-dichlore- 
thylene [17] and methyl! chlorethylene [63]. 

—Br The contribution for this group was found 
as the mean value of the increment calculated from 
methyl bromide, methyl dibromide, methyl tri- 
bromide, carbon tetrabromide [49] methyl bromo- 
acetylene [63]. 


Table 13. Sulphur containing groups 


100 Highest 


temperature 
of lit. value 


Group contribution A +- 
keal/mol 


300 600° K 600 1500° K 


A B A B 


10-942 1-026 1-167 
512 1-07 5120 107 


12-58 142 12.20 144 
14-83 051 | 
14-39 3-39 
5-58 5-26 


—I Contributions for this group were calculated 
from values of the free enthalpy functions of acetylene 
diiodide and methyl iodo-acetylene [63]. 


(e) Ring formation 

Only few data about the free enthalpy of formation 
of heterogeneous ring molecules are given in the 
literature. The only possibilities proved to be cal- 
culations of corrections for the formation of the 


0. 
and >o~ rings. The correction 


™ 


for the three-ring formation was calculated from the 
free enthalpy function of ethylene oxide [20]. 


00 ~~ 600° K 600 ~ 1500°K | of lit. value 
Group — 
A B A B K 
| | | 
HF 64.476 0-145 | 64-884 0.081 1500 
—F | —45-10 0-20 
HCL 22-100 0-215 | 22-460 0-156 1500 
—Ci | — 85 0 | — 825 1500 
HBr | — 8-723 0.234 
— Br 2-21 -0-26 
HI 0.690 0-225 
—J 15-26 0 
60 
OG? 
OL. 
l 
— 
-20 
= 
Group -60 
—SH 1500 Temp? K 
— 1500 Fig. 7. 
=s8 1500 
>so 1200 
>so, 1200 
60 
ao] | mat 
THC? 
0 
-20! 300 S00 700 900 100 1300 1500 
Fig. 6. 
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Table 15. Ring formation 


>C 


Bremner and Tuomas [9] give the heat of for- 
mation and the entropy of tetrahydrofurane at 
300° K. From this the correction for the formation 
of the 

that Uxicus rule [54] is applicable. 

The various tables for the group contributions 
show, besides these group contributions, formulae 
for the three enthalpies of formation of compounds 
differing only in one hydrogen atom from the group 
contributions. 

These formulae were calculated from data about 
CH, [39]; C,H, [39]; H,C=C=CH, [40]; H,O [39); 
H,CO[51|; HCOOH [62]; HCN [22]; 
NH, [48]; H,S [10]; HF [33]; HCI [18]; HB, [23]; 
HJ [32]. 


~ ring was calculated, assuming 


VII. EXAMPLES ON THE CALCULATION OF THE 


FREE ENTHALPY OF FORMATION 
(1) 2,2 Dimethyl propane 

C(CH,), 
300 600° K 

4—CH, 

1—C— = 

RT 


(2) 1,3 Butadiene 
H,C=C—C=CH 
2 


~43-772 
1-958 + 3-735.10%. 
0-494. 10°?. 

+ 13-089. 10°. 


2He .. = 13-960 + 2-176. 10°. 
RT ino = 0-138. 10-2. 


AG = 22-960 + 3-978. 
(3) 1 Pentyne 
HC=C—(CH,),—CH, 


—12-310 + 2436.10 7.7 
— 11-660 + 5-088. 10°? 7 
26-555 — 0-550.10°*, T 
26-700 — 0-704.10°?. T 


(4) Anthracene 


H 


600 1500° K 


RT ino 


29-285 + 6270.10 T 


2-05 +706 .10?. 
= 26-432 + 2-056 .10°?. 
0-276. 


1a 


5148 +939 .10%. 


(5) Trans 1,4 dimethyl cyclohexane 


trans 1,4 branching 


20-780 + 9-720. 
— 1410 + 5-820. 
21-886 + 4-430. 
2-80 
=— 1-128 — 1-635. 
= 1-138. 


= —48-00 + 18-47. 


(6) Tertiary butyl alcohol 
(CH,),C—OH 


600 ~ 1500° K 


3—CH, = —36-930 + 7-308.10°%. 


1—C— 


= 4385+ 3350.10 %. 


1—OH = —41-56 + 128 .10?. 


AG =—7410 + 11% 210%. 


(7) Quinoline 
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B 1—CH 
Correction A + 2—CH,— 
keal mol 
- 
Group / 
A B A B 
Som 12.86 0.63 12.86 0.63 H H 
aCy 
0. 3-53 0-16 t t 
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H H H 
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300 — 600° K 
7 HC st = 21-329 + 
= 6-608 + 0514.10. 
1 Nit =1132 +111 
AG, = 39-26 5-93 


VIII. Discussion 
Accuracy of the correlation—With the help of the 
correlation obtained in this manner free enthalpies of 
formation of gaseous compounds were calculated and 
the results compared with the free enthalpies of for- 
mation found in the literature (Table 16). It appeared 


that the values calculated in this way approach 


reality with a mean deviation of about 0-6 — , 


In the Figures] to 7 inclusive various group 
contributions have been plotted as a function of the 
temperature. Moreover—as was done in the preceding 
tables for the group contributions—those compounds 
were included in the diagrams which differ only in 
one hydrogen atom from the group contributions. 
For distinguishing from the real group contributions 
these latter materials are represented by special 
lines (~-—-—-—+—+). 

From Figure | it follows that a group contribution 
for the alkanes is, by approximation, found as half 
the sum of the preceding and the following group, in 
other words; the influence of the various hydrogen 
atoms is in these contributions approximately con- 
stant. 
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Figure 2 shows that a similar rule also holds for 
the contributions for the alkene groups, provided the 
average contributions for the H,C=C< , Ho_cH and 
po=—cH groups are considered. 

In Figure 3 the contributions for carbon-oxygen 
groups have been plotted. From this diagram the 
ready decarboxylation of the acids and esters at 
higher temperatures is apparent. 

Figure 4 shows that the stability of the nitrogen, 
sulphur and oxygen groups increases in the order as 
mentioned here. This diagram also shows the rather 
exceptional behaviour of the —O— group. 

Contributions for carbon-nitrogen groups and 
alkyne groups were represented in Figure 5. From 
this graph it appears that the entropy terms of these 
group contributions are approximately equal. 

The “‘aromatic”’ resonance groups ; 
N= and are represented in Figure 6. 

Figure 7 gives a confirmation of what might be 
expected of the —-F; —Cl; —Br and —J groups on 
account of their place in the periodical system. 

Table 17 gives a representation of the free ent- 
halpies of formation of a number of materials that 
cannot be calculated by means of this correlation and 
those of some inorganic gases. This opened the 
possibility of predicting in a first approximation the 
variation of the equilibrium with the temperature of 
the most important organic gas reactions. 


Table 17 
| Free enthalpy of formation A +- | 
100 Calculated 
Compound Formula 300 600° K 600—a°K | 

A B A B 
a C,H, 21.539 4-583 24-588 5-093 1500 39 
carbon monoxide. . . . . co 26-582 2-122 26-582 2-122 1500 39 
carbon dioxide. . . . . . co, 94-076 0-061 94-234 0-035 1500 | 39 

methanol... CH,OH 49-440 3-480 51-330 3-795 1500-63, 46, 29 
eyanogen chloride... . CNCI 34-151 0-430 34-151 0-430 1000 49 
cyanogen oxide ..... CNI 49.268 0-596 49-268 0-596 800 | 49 
ree or ee (CN), 71-889 1-059 71-889 1-059 1000 53 
sulphur dioxide ..... | SO, ~ 77-157 1-735 77-157 1.735 1500 | 16 
sulphur trioxide . ... . SO, ~ 100-203 3-943 100.203 3-943 1200 | 61 
mitzic oxide ....... NO 21-543 — 0-301 21-543 0-301 1500 16 
nitrogen dioxide... . . | NO, 7-933 1-486 7-933 1-486 1000 17 
nitrogen tetroxide . . . . | N,0, 7-112 5-460 17 
carbon oxi-sulphide. . . . | COS 39-714 0-245 39-569 ~ 0-265 1500 ll 
carbon disulphide .. . . CS, ~ 16-000 0-150 — 16-196 O177 1500 11 
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Notation 

C; D ete. — Compounds 

Elements 

K, Equilibrium constant for formation from the 
elements 

K,, Equilibrium constant 

A and B Constants 

AG, Free enthalpy of formation (American termino- 
logy: free energy of formation) 

1G? Free enthalpy of formation in the ideal gaseous 
state and at a fugacity of latm. (American 
terminology: free energy of formation) 

AG ncction ~ Free enthalpy change as a result of the reaction 
(American terminology: free energy change) 

R Gas constant 

T Absolute temperature 

Symmetry number 

? Arbitrary thermodynamic function 

al -_ Free enthalpy function in the ideal gaseous 
state at a fugacity of | atm. (American termino- 
logy: free energy function) 

H° Hg Heat content in the ideal gaseous state at a 
fugacity of | atm. 

CG Specific heat at constant pressure in the ideal 
gaseous state at a fugacity of | atm. 

AH Heat of formation in the ideal gaseous state at 


a fugacity of 1 atm. 
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Summary—The concentration of a suspension at any point is determined by measuring the capacity of a 
condenser which is designed not to interfere with the effect being studied. 

Experiments with silica powder of 56 u diameter in a mixture of benzene and nitrobenzene indicate that 
the degree of packing of the sediment is a function of its height and not of the time it is left in the settling 


tank. 


The method can also be employed for measurements in cyclone separators or in fluidized beds. 


Resumé— La mesure de la constante diélectrique d'une suspension permet de suivre l'évolution des concen- 


trations en chaque point. 


La méthode a le grand avantage de permettre des mesures instantanées, bien localisées et suffisamment 


précises. 


Si le condensateur de mesure est convenablement dessiné, sa présence ne trouble pas le phénoméne 


étudié. 


Nous avons étudié de cette fagon la décantation d'une poudre de quartz dans un liquide, dans un appa- 
reil A fonetionnement discontinu. Nos mesures indiquent notamment que le tassement du dépdét est fonction 
de la hauteur de celui-ci et non du temps de séjour dans le décanteur. 

La méme méthode de mesure pourrait s'appliquer a l'étude du fonetionnement d'un appareil de décan- 
tation continu, d'un épaisseur cyclone ou d'une colonne de fluidisation. 


Erupe DE LA D&CANTATION PAR LA MESURE DE LA 
CoNSTANTE DIELECTRIQUE 

Le but de cette note est de décrire une méthode de 

mesure rapide et précise de la concentration de suspen- 

sions de particules solides dans un liquide, a l'aide d'une 

mesure de la constante diélectrique de la suspension. 

Nous nous proposions d’étudier le probléme suivant : 
quelle est ['évolution de la concefitration d'une suspen- 
sion en ses différents points au cours d'une décantation. 
Un seul auteur, a notre connaissance, a fait des mesures 
dans ce sens (Comines [2]), il procéde A l'analyse 
d’échantillons prélevés A l'aide d'une pipette en un 
point donné de la suspension en train de décanter. 
Cette méthode présente un inconvénient: pour que 
la prise d’échantillon ne trouble pas le phénoméne de 
décantation qu'on est en train d’étudier, il faut qu'elle 
soit trés lente. Comme il ne faut pas non plus que la 
concentration varie au cours de la prise d’échantillon, 
la méthode ne s’applique qu’a l'étude de poudres dé- 
cantant trés lentement. 

La méthode par mesure de la constante diélectrique 
que nous avons utilisée donne au contraire une mesure 
instantanée et locale, ne troublant l’écoulement que 
d'une fagon négligeable. 

Le principe de la méthode est simple: on mesure 
la capacité d’un condensateur plongé dans la suspen- 
sion; le rapport de cette capacité a celle du méme con- 
densateur dans l’air donne la constante diélectrique 
de la suspension. Celle-ci est une fonction de la con- 
centration. 


En prenant un condensateur suffisamment petit et 
de forme convenable on arrive & mesurer la constante 
diélectrique en un point, ou tout au moins dans un plan 
horizontal donné du milieu étudié. L’appareil de me- 
sure que nous avons utilisé est un dieleométre (Haardt 
& Co, Diisseldorf) qui permet une mesure trés rapide 
et trés précise A la fois de capacités de faible grandeur 
(20.4 110 unF). 


D’une fagon générale la méthode présente de lin- 
térét dans tout probléme oti il s’agit de suivre l’évo- 
lution dans le temps et dans l’espace d’une concen- 
tration pouvant varier assez rapidement. 


RELATION ENTRE LA CONSTANTE DIELECTRIQUE 
ET LA CONCENTRATION DES SUSPENSIONS 


Dans les cas de suspensions de grains solides dans un 
liquide les régles de mélange simples applicables au cas 
des solutions ne sont plus valables. 


Différents auteurs ont proposé des lois dont la jus- 
tification théorique est basée sur un certain nombre 
d’hypothéses sur la forme des grains et leur répartition 
dans le milieu liquide. Ces lois ont été plus ou moins 
bien vérifiées expérimentalement dans un certain nom- 
bre de cas particuliers. 

On trouve dans l'article de BruGGEMAN [1] (1935) 
une étude théorique assez compléte des différentes for- 
mules proposées. 

&, = constante du liquide; ¢, = constante du so- 
lide; 6 = concentration de solide. 
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La formule de Rayieren [5] s’écrit : 


: 
& 
L | 
10 20 IO 40 50 
Concentration 


Fig. 1. Valeur de la constante diélectrique en fonetion de la concentration. 


La formule de LicuTENeCKER [4] est de la forme 


k dépend de la forme 


on simplifie parfois en faisant k — 0 


e= ef 


BrvuGceMaN propose 2 nouvelles for- 
mules donnant des valeurs limites cor- 
respondant aux cas de grains sphériques 


e,—e V2 


et de grains cylindriques 


2e, + & 


GUILLIEN [3] (1941) a entrepris une 
vérification expérimentale assez poussée 
de ces diverses formules, il a trouvé la 
meilleure concordance avec la formule 
de Wiener dans laquelle on fait u 
2) la formule de Licurenecker 
est souvent aussi bien vérifiée par les 
expériences pourvu que le rapport ¢, &, 


1 Loi linéaire ¢ = Oe, + (1 — d)e,. 2 Formules de Wacnnoutrz et France- 
re ne différe pas trop de lunité. 
‘o » Wac 7% et FRANCESON 
SON fy 3 Formules de ACHHOLTZ et FRANCESON (7) (1940) 
ox & tak ont également entrepris des expérionces 
en vue de contrdler les différentes 
polarisations p 4p, + 6 Formule de Rayieren formules proposées. Ces expériences 
3 sit Ge montrent que la forme des grains joue 
. 7 Formule de Wiener 
d’autant plus faible que la forme des 
Formule de Licwrenneker ¢ grains est proche de sphéres; ce résultat 
&, + 
: les conduit A rejeter les formules qui ne 
avec tiennent pas compte de ce facteur de forme; ils 
Grains Grains proposent deux valeurs limites 
spheriques cylindriques 
0 = 
2 1+a—dha + dx 
r 1 
1-85 0-306 eg (1 + — da) 
4 1+ da 
si d est petit elle devient: correspondant aux cas de particules cubiques ou 
_ Dans la formule de Wiener ils 
proposent d’écrire u = f.v avec f facteur de forme. 
La formule de Wiener [8] généralisée s'écrit : Dans le méme ordre d’idées il faut encore signaler 
— étude Voet [6] 1947 qui montre que les 
e+u & écarts par rapport la formule de BruGGEMAN pour 


ot « dépend de la forme des grains. 


des sphéres peuvent servir de mesure du degré de non 
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sphéricité des grains, et qui en déduit un moyen de 
suivre la floculation et lagglomération des grains. 

Nous avons calculé les différentes formules pour 
le cas particulier correspondant A nos expériences: 
liquide; solution de benzéne et nitrobenzene (¢— 9-94) 
solide: particules de quartz anguleuses (¢-- 4-85) les 
résultats sont inscrits sur le graphique (1) ainsi que 
les résultats des mesures directes. 

On voit qu’aucune des courbes proposées ne rend 
vraiment bien compte des nos résultats, comme le but 
de notre travail n’est pas de vérifier ces formules, mais 
bien d’étudier la décantation, nous nous sommes basés 
dans la suite sur les valeurs données par la courbe 
d’étalonnage. 


REALISATION EXPERIMENTALE 

1) Le Dielcometre—Cet appareil construit par la firme 
Haard & Co Diisseldorf est décrit en détails dans | 'artic- 
le de Wacnnoitz et Franceson. Le principe de son 
fonctionnement est le suivant: un circuit oscillant 
comporte en paralléle la capacité & mesurer et un con- 
densateur variable étalonné, on compense les varia- 
tions de la premiére en réglant le second, on vérifie 
que la somme de ces deux capacités est constante en 
comparant la fréquence du circuit oscillant a celle d'un 
autre circuit oscillant dont la fréquence est maintenue 
constante l'aide d’un quartz. La comparaison des 
deux fréquences se fait en cherchant A annuler la fré- 
quence du battement qu’elles produisent. La précision 
de la mesure dépend surtout de la précision et de l'éta- 
lonnage du condensateur variable. 

2) Le Condensateur de mesure.Aprés quelques 
essais nous avons adopté un condensateur qui com- 
prend deux armatures en forme de grille situées toutes 
deux dans un méme plan, chaque élément d'une grille 
se trouve situé entre deux éléments de l'autre. Le plan 
des grilles est orienté horizontalement, les fils, minces, 
ne troublent écoulement que d'une fagon négligeable. 
Le peu de distance (2 mm) entre les éléments successifs 
des deux armatures empéche que les lignes de force 
du champ ne divergent d'une fagon appréciable en 
dehors du plan de la grille. La mesure est done bien 
localisée dans un plan donné. 

3) Les conducteurs d'amenée— Le condensateur pro- 
prement dit étant relativement petit, sa capacité a vide 
ne dépasse pas 8 uuF la mesure est done susceptible 
d’étre fortement entachée d’erreur si on ne prend pas 
des précautions tout a fait spéciales pour éliminer les 
capacités parasites entre les conducteurs d’amenée de 
courant et les objets extérieurs raccordés A la terre. 
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Nous avons obtenu les meilleurs résultats en utilisant 
un cable coaxial dont le conducteur extérieur est rac- 
cordé A la masse et constitue blindage. Ce cAble a une 
capacité assez importante mais tout a fait invariable 
et indépendante de la position des objets extérieurs. 


ErTuDE DE LA D&CANTATION D'UNE SUSPENSION 
DE SILICE DANS UN LIQUIDE 


Le mode operatoire est le suivant—un poids donné de la 
poudre étudiée est mis dans le liquide (mélange 
opératoire de benzéne et de nitrobenzéne) contenu dans 
un evlindre gradué, le condensateur de mesure est 
fixé A un niveau donné 
et raceordé A lappareil 
de mesure, la suspension 
est rendue homogéne 
par une violente agita- 
tion du cylindre, puis 
on la laisse reposer et 
on commence immé- 
diatement A faire des 
mesures de la capacité 
toutes les cing secondes, 
jusqu’da ce que le tasse- 
ment final soit atteint. 
On recommence le 
méme processus pour 
chaque hauteur du con- 


2 Fig. 2. Le condensateur de 
reléve 


densateur, on mesure. 
ainsi des courbes sembl- 

ables a celle de la Fig.3 donnant la capacité en 
fonction du temps pour chaque hauteur, de cette famille 
de courbes on déduit ensuite facilement celles des 
Fig. 4 et 5 qui sont beaucoup plus intuitives; les cour- 
bes 4 donnent en fonction du temps les hauteurs ot 
sont réalisées différentes concentrations, les courbes 5 
donnent pour chaque valeur du temps, la concentra- 
tion en fonction de la hauteur. 


Ces courbes permettent de suivre d’une fagon quan- 
titative un processus de décantation : dés que quelques 
secondes se sont écoulées on peut reconnaitre trois 
zones: A) une zone supérieure claire dans laquelle ne 
subsiste qu'un léger voile de particules trés fines; 
B) une zone moyenne dans laquelle la concentration 
reste approximativement égale a la concentration 
initiale et C) une zone inférieure 4 concentration plus 
élevée, dans laquelle le solide est déposé. 


Alors que la simple observation visuelle permet de 
mesurer d’une fagon précise, seulement la vitesse de 
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Fig. 3. Evolution de la concentration aux différents ni- 
veaux au cours de la décantation dune poudre de quartz 
(D— 58 u) dans une solution de benzéne et de nitrobenzeé- 
ne—(Concentration initiale 21%). H initiale 20 em. 
H finale 65 cm. 
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Fig. 4. Représentation des points ou une concentration 
donnée est réalisée en fonction dui temps 
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descente du niveau A- B, la mesure de la concentration 
en chaque point a chaque instant donne une connais- 
sance compléte du phénomeéne. 

C'est ainsi que les courbes 4 mettent en évidence 
une certaine augmentation de la concentration dans 
la partie supérieure de la zone B dés les premieres se- 
condes, Cette augmentation n'est pas trés considérable 


(Lou 2%) mais elle persiste durant toute la décantation 


Fig. 5. Représentation de la concentration en fonction de 
la hauteur pour chaque valeur du temps, 


et nous avons remarquée avee des suspensions de 
concentrations diverses. C'est done un phénoméne 
assez général. 

Simultanément dailleurs la concentration de toute 
la zone B croit progressivement et lentement, cet ac- 
croissement qui atteint 2.4 3% en LOO sec pourrait 
probablement étre plus important si la décantation 
était plus lente. Nos expériences ne nous permettent 
pas de dire dés maintenant s'il en est bier ainsi. 

On peut cependant en conclure que ce n'est qu’en 
premiére approximation que cette zone B a une con- 
centration homogéne et constante. 


Un autre probléme est de connaitre l'étendue des 
zones de transition qui existent certainement entre les 
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zone A, B, C. On voit, par les courbes 4 que la con- 
centration passe de 2% A 18% sur une hauteur de 
1-5 cm et de 26% A 40% sur une hauteur de | em. Ceci 
indique done que la séparation des zones B-C est au 
moins aussi brusque que celle des zones A-B. L’ob- 
servation visuelle indique d’ailleurs que la zone de 
transition A-B n’a certainement pas 1-5 cm d’épais- 
seur. Les valeurs plus grandes déduites des mesures 
peuvent étre expliquées facilement: quelles que soient 
les précautions prises, les lignes de force du champ 
sortent un peu du plan de la grille de sorte que la 
mesure est déjd modifiée avant que la zone de tran- 
sition n’ait atteint le niveau du plan de mesure. 

De toute fagon la méthode permet de suivre facile- 
ment et avec une approximation trés suffisante la 
montée du niveau de transition B-C qu'il n'est guére 
possible d’observer directement dans la généralité 
des cas. 

Les mesures de constante diélectrique nous per- 
mettent encore de mieux connaitre le probléme assez 
controversé du processus de tassement de la zone C. 

On constate en effet que le dépét n’atteint pas 
demblée la concentration maximum qu’il est susceptib- 
le d’atteindre finalement. Certains auteurs sont d’avis 
que le tassement maximum n’est atteint qu’au bout 
d’un délai relativement long au cours duquel les parti- 
cules glissent l'une sur l'autre, le mouvement étant 
entravé par les frottement entre grains. Pour ces 
auteurs la concentration de la boue épaisse obtenue 
dans un décanteur industriel dépend du temps de 
séjour dans le décanteur. 

D’autres auteurs pensent que le tassement est da 
a la compression exercée par le poids des couches de 
dépét successives. A chaque pression correspond un 
état d’équilibre entre les grains qui est rapidement 
atteint. Pour ces auteurs la concentration de la boue 
d'un décanteur dépend de la hauteur de la couche de 
dépédt. 

Les mesures que nous avons faites montrent que 
le tassement continue 4 augmenter en un point de la 
zone C tant que celle-ci augmente c’est-a-dire tant que 
le poids de poudre déposé augmente, mais ce tassement 
cesse de se modifier dés que ce poids reste constant. 

Ces mesures montrent aussi que le tassement final 
varie d’une fagon appréciable d’un niveau A l'autre de 
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la zone C, en diminuant de bas en haut, ceci confirme 
que le facteur principal dont dépend le tassement est 
la compression due au poids des couches supérieures 
le temps nécessaire pour réaliser |’équilibre pour un 
état de force donné étant inappréciable. 

Il faut remarquer ici que les expériences ont été 
faites avec une poudre relativement grosse (silice 
C:dm = 58u) ne floculant pas il est probable que le 
facteur temps est plus important quand il y a flocu- 
lation. Cette hypothése explique d’ailleurs le désac- 
cord entre les différents auteurs. 


Finalement il faut encore indiquer que la méthode 
de mesure est surtout intéressante pour suivre |’évo- 
lution des concentrations plut6t que pour en faire des 
mesures absolues trés précises. I] est en effet toujours 
possible que la présence de la grille modifie légérement 
les concentrations, cela est méme probable lors des me- 
sures dans la zone C plus ou moins tassée, mais il est 
bien certain que le sens et la durée des variations 
de ces concentrations n'est pas changé. 


Le présent travail a été effectué au laboratoire de 
Chimie Industrielle de l'Université de Bruxelles sous 
la direction de Monsieur le Professeur pe KrysEr 
grace A une subvention du Comité Tassel. 


NOTATION 


e — constante diélectrique de la suspension 
é, — constante diélectrique du solide 
&, — constante diélectrique du liquide 
6, — fraction du volume occupé par le solide 
6, —— fraction du volume occupé par le liquide 
d,, — diamétre équivalent moyen des grains 
H, ~ hauteur initiale de la suspension 
H, — hauteur finale du dépét de poudre. 
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Calculation of the efficiency of counter-current stage-wise mass transfer processes 
with constant distribution factor, when in the stationary state * 


I. Distribution of one component only 
A. KLINKEN BERG 


N.V. De Bataafsche Petroleum Maatschappij, Royal-Dutch/Shell Group, The Hague 
( Received 19 July 1951) 


Summary—A rapid and easy method is presented for the derivation of equations for counter-current 
stage-wise mass transfer of a single solute with a constant distribution factor, when in the stationary state. 
The method is valid for a group of related processes but it is specifically applied to extraction. 

Equations are derived for extraction with an intermediate feed point (“complete rectification’’) and for 
extraction with introduction of the feed at an end (“‘extraction” or “ washing’). Some of their consequences 
are discussed. 

The use of reflux is included in the treatment. 

Equations are also given for the accumulation of solute in the feed stage. 

The separation of two solutes is to be discussed in a subsequent paper. 


Résumé—L ‘auteur présente une méthode rapide et commode pour la dérivation d'équations convenant 
& l'état stationnaire dans les problémes de transfert de masse étagés A contre-courant, lorsque le systéme 
présente un facteur de distribution constant. Cette méthode, valide pour tout un groupe de procédés 


d'échange, est particuliérement applicable & extraction. 

Dans ce cas, les equations s'appliquent aussi bien & la rectification compléte (alimentation é un niveau 
intermédiaire) qu’é l'extraction simple ou lavage (alimentation en téte). 

L'auteur discute certaines conséquences de ces équations. I] étudie également le cas du reflux. I] donne 
aussi les équations correspondant l'accumulation de corps dissous dans I'étage d'alimentation. 

La séparation de deux corps dissous fera l'objet d'une publication ultérieure. 


1. Diusrrisution Facrors as THE EXTRACTION 
Facror 


The following calculations apply to counter-current 
mass transfer processes, carried out in ideal stages 
under such conditions that the amounts of a certain 
solute leaving a stage with the two streams have a 
constant ratio. 

In counter-current extraction the amount of the 
solute, whose extraction is being studied, may be small 
with regard to the amount of liquid originally carrying 
this solute and also with regard to the amount of 
solvent. Under these conditions we may neglect the 
volume effect of the extraction so that the throughputs 
of the two streams will not change from stage to stage. 

Under these same conditions in most cases the ratio 
of the concentrations of the solute in the two streams 
will also be constant (NeRNsts Distribution Law). This 


For the purpose of this paper throughputs are much 
more important than concentrations. 

If the concentrations (C’, ¢c) are multiplied by the 
throughputs of the phases (volume per time): 

V = throughput of solvent rich or extract phase 
in volumes per unit time; 

v = thronghput of solvent poor or residue phase in 
volumes per unit time; 
we obtain 

Q = throughput of solute by solvent rich phase 
(any unit per unit time); 

q = throughput of solute by solvent poor phase 
(same unit per unit time). 

Under equilibrium conditions, i.e. when leaving an 
ideal stage, Q and gq will also have a constant ratio, the 
extraction factor E. 


_ amount of solute leaving ideal stage with extract phase 


constant ratio is called distribution coefficient or “amount of solute leaving ideal stage with residue phase ‘~) 
partition coefficient K. The concentration is to be. ¢hat 
expressed in any units per unit of volume. K (3) 


We define 
K € __ concentration in solvent rich phase 
¢ concentration in solvent poor phase 


(1) 
under equilibrium conditions. 
* This paper was presented to the XII International 


Congress of Pure and Applied Chemistry, New York, September 
10-13, 1951. 


q ev v 


The same reasoning applies to the class of extraction 
provesses sometimes called “fractional distribution,” 
where a mixture of solutes is added to two solvents 
flowing in counter-current. However, the choice of the 
phases to be used in numerator and denominator in 
the above definitions becomes an arbitrary one. 
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The concept extraction factor may be generalized to 
cover other unit operations. We may call it in general 
a distribution factor. Such distribution factors are 
of considerable use in the study of: 

Liquid-liquid extraction. 

Solid-liquid extraction (leaching). 

Stripping. 

Absorption (under isothermal conditions), 
where they are often very nearly constant. 

In distillation the distribution factors in general 
are variable because of the temperature gradient. 
Fenskes well-known equation, however, is of the type 
discussed here. 

It should be realized that the low concentrations re- 
quired for the application of constant distribution fac- 
tors should not occur in extraction with two solvents. 
If they do, the system has been insufficiently loaded 
with solutes and is, therefore, not working econom- 
ically. 

It is felt that calculations based on constant dis- 
tribution factors will: 

1. in many cases offer a good quantitative approx- 
imation ; 

2. very often be useful for rough estimates ; 

3. be extremely useful for educational purposes, 
because many features of the stage-wise mass transfer 
processes can be studied here in their simplest form. 

The perhaps oldest application of distribution fac- 
tors is found in absorption, where Kremser [1] in- 
troduced the absorption factor A. 


2. DEFINITION OF PROCESSES TO BE STUDIED 


The processes to be considered are continuous in the 
time. We shall not consider working on a batch of feed, 
as in batch rectification, nor discuss transient pheno- 
mena occurring before the stationary state is reached. 

They are operated in stages, i.e. the concentrations, 
etc. assume a number of discrete values as a function 
of another discrete variable, the stage number. 

The stages may be operated either intermittently, 
as in separating funnels, or non-intermittently, as in 
systems of mixers and settlers. This difference does 
not influence the results of the calculations. 

The counter-current principle as such does not re- 
quire clarification. However, it is to be stated parti- 
cularly that the feed need not be introduced at an end 
of the system. A treatment will be developed which is 
valid for two solvents flowing counter-currently with 
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introduction of feed mixture at some intermediate 
point. 

In continuous distillation it is obvious that the 
mixture to be distilled must not be introduced at the 
top (or the bottom) of the column, since this would 
not allow obtaining the more (or less) volatile compo- 
nent in the pure state. Hence the use of columns for 
complete rectification, comprising a “rectifying sec- 
tion” and a “stripping section.” 

The same object was achieved in extracvion by VAN 
Disexk [2] by the use of two imn.iscible solvents, flow- 
ing in counter-current with introduction of the feed at 
some intermediate point. This process is in use for 
the extraction of lubricating oils (‘Duo Sol Process,” 
solvents: propane and a mixture of phenol and cresylic 
acid) and citrus oils (vaN Disck and Rvuys [3], van 
Wuk and van Dick [4], solvents: pentane and aque- 
ous methanol). More recently Scuerpet [5] called 
attention to it. Here again we shall speak of complete 
rectification. Rectification is understood to be a pro- 
cess for purifying a phase by counter-current contact 
with another phase containing certain materials (the 
reflux) produced from the first phase. Obviously this 
definition applies on both sides of the feed stage, in 
distillation, extraction and several other processes. 

In distillation the reflux is usually produced out- 
side the system (in condenser and reboiler) ; in extrac- 
tion it is usually—when both solvents are originally 
free of solute—produced within the system. 


3. Basic Equations 
The calculations will employ extraction factors and 
throughputs of solute. 

The following convention is used for defining ex- 
traction factors with the aid of eq. (2) when a polar 
solvent is used in counter-current to a hydrocarbon 
solvent: 

phase rich in polar solvent is called extract phase, 

phase rich in hydrocarbon solvent is called residue 
phase. 

The section of a complete rectification between feed 
and inlet of polar solvent by analogy of single-solvent 
extraction is called an extraction section; the section 
between feed and inlet of less polar solvent a washing 
section. 

This is shown in Fig. 1. 

The stages will be numbered in the direction of the 
flow of the non-polar solvent; the transport of solute 
thereby is represented by the q values. 
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Flow rates of solute carry a subscript to show the — with the aid of (4) is rewritten in the form 
stage of origin. 
A section of a counter-current system is represented wlth 
by Fig. 2. The squares: represent the stages*, the oF, after multiplication by — 
arrows the flows of solute. Q, + + 0. (8) 


the extraction in the &'th 


This equation Jescribing 
stage is found by giving the 


streams entering the stage 
section section a positive and those leaving 

extract phase | polar solvent negative sign, multiply- 
tec? ae = ing each stream by a power 
Fig. 1. Definitions of various terms. of EB, the exponent of which 

is the number of the stage 

For a system of ideal stages in stationary operation into which the stream enters, and equating the sum 
the basic equations are: to zero, 

Upon adding the corresponding equation for the 
adjacent stage (k + 1) the terms for the transports 
de = Constant (material balance). (5) between stages k and k +1 are seen to cancel. This 
procedure may be continued in the same way. 


Eq (equilibrium condition), (4) 


The latter equation is a finite difference equation. ; 
Consequently for any counter-current system with 


constant extraction factor the procedure for finding 
the equations governing the extraction is the following: 


This set of equations may be solved by starting 
from the end(s) of the system and working towards 


Skeet 9 1. Number the stages in the direction of the flow 
| resides phase of that phase whose solute flow is in the definition of 
| in the denominator (i.¢. the residue phase). 
2. Assign a positive sign to solute flow rates enter- 
wal T extract phase ing, and a negative sign to flows leaving the extrac- 
QO, tion system. 
Fig. 2. Definitions of flow rates. 3. Multiply each solute flow rate by a power of E 


the exponent of which is the number of the stage in 
the feed such as was Kremsers procedure when he which such flow enters. 


derived an equation describing counter-current ab- For this purpose it will be necessary to add ficti- 
sorption [1]. Scnetser (5) when studying the process tious stages for all streams leaving the system. 
of complete rectification in extraction proceeded 4. Equate the sum of the flows found under 2 
numerically from both ends toward the feed. to zero (= overall material balance). 

The calculus of finite differences has been employed 5. Equate the sum of the products found under 3 
by and Tour [6] and [7], who showed zero (= extraction equation). 
that the solutions assume the form of series of powers 6. Reduce the two equations to a single one by 
of the extraction factor. eliminating one of the unknowns. 


The present paper tends to avoid the lack of ele- 
gance of the step-wise procedure on the one hand and 
the use of a not too familiar mathematical technique 
on the other. 


To this end the material balance of the é'th stage 


4. AppLication To CompLere RecriricaTion 
Fig. 3 represents a complete rectification, inn —m —1 
stages, with feed F in the m’th stage. 

According to the procedure described in section 3 
the two equations are found to be: 


In a conventional drawing of a plate column the divid- - 


ing lines represent stages. —Q, + + F + Qaim = 9. (10) 
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Case 1.—No reflux at the ends 


If the solvents are pure, so that no reflux is applied 
at the ends, we have: 


(11) 
Crim = 9, (12) 
so that: 
1 
=F (14) 


and the ratio R of the amounts of solute in the two 
products is given by: 


(E" —1) 
(15) 
Special cases are 
a) E=1 
m Yo 
(16) 1 
R= (18) 
b) 
Rx k* (19) 
E< 
(20) 
qen-1* F +1 (21) 
R= E*. (23) 


The symmetrical case discussed under d has been 
studied by Romersen [8]. 

Equation [15] for the general case has been derived 
by Barres and Kieman [9]. 

Related equations have also been derived by STENE 
[10], but his mathematical approach, through inter- 
mittent batch processes (non-stationary extractions in 
separating funnels), is rather complicated. 


Case 11.— Total reflux 

Total reflux is applied when the solute leaving the 
system at either end is re-introduced with the other 
solvent at that same end, i.e. 


Yo = 
Qaim = 
Equation (9) then gives: 


qnim-i 
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In a symmetrical system (n — m) the ratio R 
obtainable is: 

E*"~' with total reflux (eq. 24); 

E" without external reflux (eq. 23). 

As Romertscu [8] has pointed out, this shows what 
maximum improvement is obtainable by the use of 
reflux. 


Case 111 .--Finite reflux 


Reflux ratios r, and r, at the ends of the extraction 
and the washing section are defined by: 


qusm-1 
Qn le 
re +1 
fo 
F 
| 
| 9n+m-! 
| m nem-! 
| <—— nem 
nstages 
Fig. 3. Complete rectification. 
Substitution in eq. (9) and (10) gives: 
R— qo 
qnim-1 Qn +m (25) 
(re + 1) 


_ 
(re + 1) E™ + re EB — (te + 1) | 

Use of reflux has been discussed by AsseLin and 
Comincs [11], who have used graphical methods of com- 
putation. They stated that it may improve the separa- 
tion obtainable; however, the higher solute concen- 
trations will counteract this through a reduction of the 
separating effect of the single stage. The accumulation 
of solute which is responsible for this reduction will 
be discussed in Section 8. The use of reflux in the 
washing section has been described by Tismstra [12]. 


Discussion of Case I, no reflux 
The significance of the equations (21), (22) and (23) 
for symmetrical systems is easy to visualize. Eq. (23) 
is represented in Fig. 4. 

Thus, if Z = 2 in a system of II stages, n = 6, 
the ratio R = E" = 64 so that 985% of the feed 
leaves the system with the extract phase and 1-5% 
with the residue phase. 

In a plot of log R against log EZ a straight line is 
obtained with an angle of inclination «, for which 
tga =n. 
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Such a graph is shown in Fig. 5. 5. Appiication TO ORDINARY COUNTER-CURRENT 
EXTRACTION WITH AN INITIALLY PURE SOLVENT 


Asymmetrical systems, described by eq. (15), are 
represented in this graph by curves having asymptotes Equations for extraction with a pure solvent in an 
under inclination tg x = _m [see eq. (20)] and tg —» extraction section only (see Fig. 6) are obtained by 
[see eq. (19)}. substituting m <1 in eq. (13)-(15), or by straight 
application of the principles of section 3. 


E(E"—-1) 
PF i’ (27) 


E* 


+ 


E(E"-1) 


1 


known as equation. 


Equation (26) is 

The consequences have been discussed by many 
authors and it is found in various other forms, see ¢.g. 
Perry {13}. 


It is readily seen that the extraction can be made 


0 | | | | to approach completeness, Lim gq, 0, only for 
0102 05 1 2 Ss 10 Eel. on 
E 


Fig. 4. Fraction extracted in a system for complete recti- For E11, Lim dn * Fi(l—E£). In words: If 


fication consisting of 2 n-1 stages with equal extraction ‘ 
factors and in which the feed is introduced in the middle. the solvent is not good enough (K too low) or not 


enough of it is used (l/r too 


This is also easily seen from 
the graph Fig. 7. 


KReEMSERs equation allows 
solving such problems as deter- 
++ mining the minimum amount of 
solvent required (corresponding 
to E = 1) and determining the 
economic optimum for the com- 
it bination of cost of solvent cir- 
HH culation and cost of stages. 


10} +4444 


6. APPLICATION 
TO A WasHIne Process 


The washing process (see Fig. 8) 
is obtained by putting = 1, 
whence from eq. (13)-(15): 


a01 10 100 


Fig. 5. Ratio in which the feed is divided over the two products in_ various R— EV (El) 
systems with constant extraction factor. 


(31) 
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7. Appiication TO OrpiInaRY CoUNTER-CURKENT 
EXTRACTION WITH AN ImMpURE SOLVENT 
When the solvent is impure (Q,., +0) the two 

fundamental equations read: 


+P 
—Q, + FE +Q,.,2"—q, =0 


whence 
| 


In =F l + Betta (32) 


q» is the sum of two terms representing 
the result of an n stage extraction of the 
feed F with pure solvent [eq. (26)] and an 
n stage washing process carried out on the im- 
purity [eq. (29)]. 

Equation (32) may be re-arranged to: 
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E(E"-1) 


This ratio approaches unity for E > 1 or E< 1 and reaches 
a maximum value equal to n for B= 1. 

Accumulation, therefore, is most pronounced for the com- 
ponent with £ = 1 and in systems with many stages. 

The accumulations are of practical importance because 
they increase the miscibility of the solvents, thereby reducing 
the selectivity of the separation, so that they may limit the 
throughput. 


| 


Fig. 6. Extraction section only. 


mea ner) & as 


Hence the alternative method of describing 
this process consists in deducing the solute 
flow in the residue phase that would be in 


equilibrium with the solute flow @,., in the 


extract phase (viz. 


subtracting this from 

all solute flows in the residue phase and apply - 

ing KREMSERs equation to the remainder. 
Both forms of the equations have been given 


by Soupers and Brown [14]. 


+ 


= 


0 | Li 
alos 05 1 2 #5 10» 


Fig. 7. Fraction extracted in an extraction section consisting of 


8. ACCUMULATION OF SOLUTE 


n-stages with equal extraction factors. 


The feed stage (m) in a system for com- F 

plete rectification (Fig. 3) receives flows | 

and Q,,., from adjacent stages, so 

that the solute load may be a multiple 
of that caused by the feed only. o<— aici 


The total load of the feed stage is 
also given by dm + Qm == (E+ 1) dm - 

One may calculate q,, from and Q,. » 
by applying equations for extraction with impure 
solvent. 


The following is found: 
qm + Qu E+1 x 
F 


+ re (Tw + 1)} \(re + ve (34) 


Conclusions worth noting are: 
1. For a symmetrical system (nm) in the absence of 
reflux (r, = 0; 0): 


qn + Qn E+1 k*—1 


(35) 


~ 


Fig. 8. Washing section only. 


They also increase the time required for reaching stationary 
conditions. 

2. If BE” and E” > 1, when the ratio R of solute in final 
extract phase to solute in final residue phase according to 
eq. (25) is very high, eq. (34) reduces to: 

qm + Qm = E+1 
F E—1 
ie. the accumulation is independent of the reflux ratios. 
Similarly for low R, <1; <1 
qm + Qm E+1 


F 


Consequently if a sharp separation is obtained between com- 
ponents | and 2 with 2, < 1 and R,»> 1 the separation can 


(35a) 


(35 b) 
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be improved by the use of reflux, whereas the accumulation 
of solute and therefore the maximum throughput are, in first 
approximation, not changed. 

The accumulation has also been studied by STene 
[10] for systems of separating tunnels. Equation (35a) 
due to Stene, has been used by O'Kerrre, DoLLiver 
and Sritier [15] for the study of the accumulation 
of two streptomycins and “B") in the feed 
stage of a symmetrical system. In their experiment 
the distribution coefficients were Ay = 2-2; Ay 
so that Ay Ay 
justed so as fo make Ey 
1 | 3-49. 

The ratio 
was found to be for both components 2-15. Equation 
(35a) also gives 2-15. From eq. (35), the exact equa- 


3-49. The ratio Vie was ad- 
| 349 


I-87 and Ky, 


in all 11 stage system (m — 6) 
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tion to be applied in this case, a value 2-05 is found. 
Obviously for E = 1 Srenrs equation, as quoted by 
O'Keerre et al. would not be applicable at all, Stenr, 
however, also gives the exact eq. (35). 

The rate of approach to the steady state in sym- 
metrical systems has been discussed by SCHEIBEL [6]. 
He showed that the approach to steady conditions is 
slower the higher the number of stages, and further- 
more, it is slowest for the component with 2 = 1. 
This is understandable in view of the similar trend of 
the accumulation. 


9. TO THE PROBLEM 

OF THE SEPARATION OF TWO SOLUTES 
Two solutes | and 2 are characterized by their dis: ri- 
bution coeficients A, and K,, where K,/K, =p 


(8 > 1). 


NOTATION 


Definition 


concentration in solvent-rich phase 
concentration in solvent-poor phase 
extraction factor [see eq. (2)] 
rate of input of solute as feed 
distribution coefficient [see eq. (1)]} 


number of stages in washing section (incl. feed stage) 
number of stages in extraction section (incl. feed stage) 
rate of throughout of solute by solvent-rich phase 


rate of throughput of solute by solvent-poor phase 


ratio of the amounts of solute in final extract and extraction 


residue [see eq. (15)] 
reflux ratio at the end of the extraction zone 
reflux ratio at the end of the washing zone 
throughput of solvent-rich phase 
throughput of solvent-poor phase 
angle of inchination in Fig. 5 
ratio of distribution coefficients K, A, 


Unit 


| any unit per unit volume 
any unit per unit volume 
dimensionless 
same units as 
dimensionless 
dimensionless 
dimensionless 
same unit as used in defining concen- 
trations, per unit time 
same unit as used in defining concen- 
trations, per unit time 
dimensionless 


dimensionless 

dimensionless 
volume per unit time 
volume per unit time 
dimensionless 
dimensionless 


Indices indicating stages (see Fig. | and 2): 


fictitious stage beyond end of washing section 
last stage of washing section 

any stage 

feed stage 

last stage of extraction section 


fictitious stage beyond end of extraction section 


Indices indicating components: 


component with lower distribution coefficient 
component with higher distribution coefficient 
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In view of eq. (3) their extraction factors will have 
this same ratio: 


The factor § characterizes the difficulty of the 
separation. It is equivalent to the relative volatility 
in distillation. 

We may now be faced with either of two problems. 

1. Given: Feed composition, partition coefficients 
A, and Ay, conditions of extraction (number of stages, 
feed point, and ratio of phase volumes). 

Asked : 

Since n,m, Ey and Ey are known, this is easily 
solved with the aid of eq. (13) and (14). 


Product compositions, 


2. Given: Feed composition, product compositions 


and ratio § of partition coefficients. 
Asked: Conditions of extraction (mn, m, By, Bs). 
The three given compositions may also be replaced 


by two compositions and a vield, 


1951 \. H. A, Lacwrnrer and 
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This second problem is of considerably more prac- 
tical importance and is also more difficult to handle. 


It will be exammed in the second part of this paper. 
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Calculation of the efficiency of counter-current stage-wise mass transfer processes 
with constant distribution factor, when in the stationary state + 


Il. Stage requirements for the separation of two components by two solvent extraction 


A. KLINKENBERG *, 


H. A. Lauwerter ** and G. H. Reman ** 
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Summary 


The equations governing counter-current extraction as derived in Part 1, “Distribution of one 


component only,” are applied to the two components of a binary mixture. 


Methods are described to derive the minimum number of stages required for the separation of a given 


feed into products of a given composition, 


The results were applied to experiments reported by ScHEIBEL. 


Resumé 


Les equations correspondant & Uextraction & contre-courant, établies dans la lére partie, sont 


appliquees aux deux constituants d'un melange binaire. Les auteurs décrivent la méthode permettant de 


determiner le nombre minimum d'étages nécessaires pour la séparation d'une alimentation donnée en deux 


produits de composition imposée. Les résultats sont appliqués & des résultats expérimentaux de SCHEIBEL. 


A problem generally encountered in practice is that 
of separating a given feed by two-solvent extraction 
into two products of specified composition while the 
coefficients of distribution of the two components for 
the solvents used are known. The question is to find 
the total number of stages, the feed point and the 
extraction factors of the two components. 

In the literature special cases of this problem have 
been treated on several occasions. 

Van Disck and Scuaarsma [1] obtain a good 
separation of a mixture of equal amounts of pro- 


acid and propane-1,3-dicarbonic 
acid by means of extraction with ethyl acetate and 
water as solvents in a symmetrical extraction system 
under conditions that the product of the extraction 
factors of the two components equals 1. This result 
STENE [2] for 
In a symmetrical 


agrees with calculations made by 
symmetrical extraction systems. 


+ This paper was presented to the XII International 
Congress of Pure and Applied Chemistry, New York, September 
10-13, 1951. 

* N.V. De Bataafsche Petroleum Maatschappij, Roya 
Dutch Shell Group, The Hague. 

** Koninklijke/Shell-Laboratorium, Amsterdam. 
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system a mixture of equal amounts of two components 
is split into products of equal purity when the pro- 
duct of the excraction factors equals 1. 

Busia and Densen [3] treat the problem of the 
separation of a mixture of different amounts of two 
components according to Crates procedure. They too 
mention that for a symmetrical system and for pro- 
ducts of equal purity the product of the extraction 
factors of the two components should be equal to 1. 

The treatment of Romerscu [4] is likewise restricted 
to symmetrical systems, 

However. no information is found in the literature 
on the general problem as stated above. In the follow- 
ing the mathematical formulation of the problem and 
a numerical as well as a graphical method for solving 
the equations will be given. 


2. Discussion OF SPECIAL SOLUTIONS 

For the treatment of the problem we start from the 
ratios R, and R, in which the two components are 
distributed over the two products. These quantities 
are readily calculated from the overall material 
balances for both components. In this respect it is 
sufficient that either the compositions of the feed and 
the products are known or two of these three compo- 
sitions together with the vield, 

As explained in the preceding paper* the following 


equations hold good for R, and R,: 


Ew 1) 


R,, 
(EY 
R,. (2) 
Further: 
KE. p (3) 
“1 


Component 2 is to be the one with the higher 
distribution coefficient, so that 8 > 1. Then in practical 
Ry» 1 and 1. In order to 
measure the E-values it 


cases E,>1 
is further convenient to 
introduce the parameter 

log By 


log (Sa) 


(see graphical method). 
As a consequence of the above: O<2< 1. 

In the three equations given R,, R, and f are 
known, it is asked to find Vv, hence the E-values, 
n and m. These 3 equations with 4 unknowns still 


* Same title, Part 1, Distribution of one component only, 
by A. Kurskenserc. Chem. Eng. Science 1951 1 86. 
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allow an infinite number of solutions. It is reasonable 


to demand that » = m— 1, the total number of stages, 
be a minimum. 


These simultaneous equations of higher, and possib- 
lv fractional, order do not allow a. straightforward 
solution. In several special cases, however, the equa- 
tions can be simplified, 

These include : 

Casel x 

p'. (da) 

(4b) 


Method of solution: 


From (2) and (4b) 


(9) 


From (1), (4a) and (5): 


(3 Rom ) 


R th 


Find m by trial and error, and then « from m 


and (5). 


( ‘ase 


2m " (7) 
Method of solution: 
From (1) and (7): 
BP 
From (2) and (7): 


From (3), (8) and (9): 


| 

Solve for m and then find » from (7). 
from (8) and then from (3). Calculate xr. 


(10) 
Find E, 
Case 3. x 


(lla) 


ky 

Method of solution: 
Divide (2) by (1) and substitute (Ila) and (ILb): 


R, 


VOL, 
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Substitute (12) in (2): 


R, 14+R 
Sm 2 1 
tars % 
1+R 
» 2 
m 
log 
Find » from (13) and (15): 
2 log | 
log B 
Case 4. 
m 


Method of solution: 


From (1) and (2) after substituting m — n: 


EY = R,, 
R,. 
Dividing (19) by (18): 
R 
OL. 
952 
log R, 
"log 
k, 
2 


(14) 


(15) 


(16) 


(17) 


(21) 


(22 


Find £, and EB, from (18), (19) and (21). Calculate x 


Case 5. 
m 2n. 

Method of solution: 
From (1) and (23): 

R,. 
From (2) and (23): 

Ey Ry + Ry. 
From (3), (24) and (25): 


Solve for » and then find m from (23). 
from (24) and then F, from (3). Calculate x. 


Case 6. x l 


Method of solution: 
From (1) and (27a): 


23) 


24) 


(25) 


(26) 


Find 


(27a) 


(27b) 


(28) 
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From (2), (27 b) and (28): 


pm 


Find m by trial and error, and then # from m 


1) 


7% (29) 


and (5). 

It is to be noted that cases 3 and 4 require the same 
total number of stages, so that the optimum number 
corresponds to conditions intermediate between those 
of these two cases. 

A plot can be made showing (n — m— 1) against 

log By 
log 
The above theory has been applied to experimental 


results obtained by ScneiBer [5| for the separation 
of ortho- and parachloronitrobenzene with a heptane 
mixture (“Skellysolve C*\) and aqueous methanol 
(15% water) as solvents. 

There was a slight change (of 2°) in the extraction 
factors at the feed inlet. This was neglected in the 
present treatment. 

For this pair of substances in the solvents men- 
tioned EB, = ortho 
and component | 


1-620, where component 2 
para and 


E amount of component in methanol ; 
amount of component in heptane 

In expt. 5 the product quantities and compositions 
were: in methanol solution 3-39 lb/hr with 12% para 
isomer (solvent-free), in heptane solution 1-96 Ib hr 
with 81-7% para isomer, so that: 


256 
R, (para) 1-96 
3-30. OSS 
“28. 
R, (ortho) 1-96 


The results obtained are given in the following 
table, together with the #-values corresponding to 
NCHEIBELs experimental conditions and the stage 
numbers calculated by his trial and error method. 
The stage numbers have been plotted against « in 
Fig. 1. 

It is seen that the total number of stages (” + 
m — 1) passes through a flat minimum. Consequently 
the total number of stages corresponding to both the 
symmetrical position of the feed (case 4) and the 
symmetrical position of the #-values (case 3) will in 
many cases be sufficient as an estimate of the minimum 
total number of ideal stages. 

It is also seen, upon comparing cases 3 and 4, that 
a change in position of the feed under nearly optimum 


A. H. A. Lacwerten and G. H. Reman: Calculation of the efficiency. I 


Table 1 


E, 
(para) (ortho) 
0 baw 
Od 0790 1-27 
134 
1-00 1-62 


SCHEIBEL 


expt. 5 0-627 0836 1,353 


(mean) 


SCHEIBEL EXP $ 


5 


10 
og 


Fig. 1. Evaluation of Senetmens expt. 5. 


ae 


conditions can simply be corrected by a slight change 
in the solvent ratio. 

Safe operating conditions will be in the range be- 
tween cases 3 (x = 4) and 4 (n = m). Outside this 
range the total number of stages increases rapidly. 

SCHEIBELs expt. 5 corresponded to conditions 
only just outside the range mentioned. 

SCHEIBELs total number of ideal stages was 13-6. 
Cases 3 and 4 require 13-4 ideal stages and it is seen 
from Fig. | that the minimum could not be below 
about 13-2. 

Further it is to be noted that for cases 3 and 4 the 
total number of stages depends on the ratio of R, to 
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R, only, (.e. only product compositions are required 
and feed composition and yields have no influence. 

The above method does not allow of an accurate 
caleulation of the minimum number, and the value 
of m-- n— 1 obtained for cases 3 and 4 (viz. eq. 22) 
would have to be used as an upper estimate. 

A graphical method was therefore worked out which 
allows finding all solutions of the original equations 


instead of a limited number only. 


3. A GENERAL GRAPHICAL SOLUTION 

The following equations are given [see eq. (1)-(3)]: 
— + R,) 4 

— ET (1 + R,) + 


(30) 


They can be brought in the following symmetrical 
form 
Ey "a R,) 


Ey R,) 


+ 
Introducing new variables by means of 
| 
nv) log (m4 n) log B | 


we find 


(31) 


R, 
Ry 


or finally, when instead of R, and R, new parameters 
a and b defined by 

1+ Rk, 

are used, 
wre 
(1 
r+ye=l. 


a 


(33) 


From (32) it follows that a and 4 are in general small 
quantities. The physical meaning of a and + is clear; 
they denote the small amounts of the one component, 
relative to the feed, in the final stage, where a large 
amount of the other product is obtained. 

The minimum problem determined by (33) and 
* — minimum can be solved graphically as described 
below. 

Since the system (33) is independent of £ the 
following important theorems are derived from eq. (31). 

a) The total number of equilibrium stages plus one 
required is inversely proportional to the logarithm of 
the ratio # of the distribution coefficients. 


VOL, 
l 


Case m ” nam 
34 304 
2 | lel 42 
3 62 x2 3-4 
4 7-2 7-2 3-4 
6 Is 405 21-5 
NUMBER OF STAGES 
gE, 
+ + + + 4 + 4 1952 
+ + o + + 
feom.? 
% 
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b) The location of the feed as determined by the 


ratio is independent of 


c) The logarithms of the extraction factors are 


proportional to log f. 

The graphical solution proceeds as follows 
(Fig. 2): 

In a (w, t) diagram the following two curves 
are drawn 


C,: w log {(1 — a) 10-6 + 108 
Cy w log {(1— b) 4 b 108. 


ivery line / through the origin determines a 
solution of (33). Let / meet C, and C, in the points 
P(t) and P,(t,); then we obtain, if the slope of / 
is represented by @; from equations (33): 


The middle M of P, P, has the abscissa 5s. 
Consequently, the extreme left-hand point M, of 
the curve C, described by + 
M for variable 0 deter- 
mines the minimum total 
number of stages. When 
the values of s, O, t, and 
ty corresponding to the min- 
imum are denoted by «,, 
t and we get the 
complete solution in the 
following form 


For a rather wide range 


of the values of a and + 
the C,-curve is found to be 


very flat in the vicinity F 


of M,. In the case of an 

error in @, (i.e. the location of the feed point) the 
value of & is not appreciably influenced. However, 
f} and & (i.e. the extraction factors) are very suscept- 
ible to a slight alteration of O,. Obviously, an error 
in O, has no effect on the purity of the end products, 
the purity being determined by the choice of the C, 
and curves. 


These conclusions of course, are in line with those 
found in section 2. 


To facilitate calculations, a graph (Fig. 3) has been 
constructed containing a set of C, and Cy curves for 


Fig. 2. The graphical determination of the total number of the 
stages, the feed location and the solven tratio with the aid of a 
w, ¢ diagram. 


CURVES FOR VARIOUS 0 VALUES 
CURVES FOR VARIOUS VALUES 


ig. 3. Nomogram of C, and C, curves. 


different a and 6 values. As this graph is symmetric 
with respect to the ¢ axis, only the upper part is given. 
Hence, the dashed C, curves appear as the reflected 
continuations of the corresponding C, curves. For any 
special extraction problem the corresponding C, and C, 
curves have to be selected and the part of the C, curve 
‘in the vicinity of the minimum position has to be 


SA Bm 
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constructed 
lines 
required information can be obtained. 


A set of curves (Fig. 4) has been constructed giving 


the values of s, for different a en > values. 


n) log curves in an (a,b) 


Lim 4 


Fig. 4. Graph of «, 
diagram Chart of relative minimum total number of 
stages. 


Thus if 80% of component | leaves with the final 
residue phase and 90% of component 2 with the extract 


phase, a = 0-20 and 4 = 0-10, so that s = 1-55. 

“ 
92 

e-090%8 6202038 

ot 

02 03 os 06 07 06 os 10 


Fig. 5. Determination of minimum number of stages for product 


compositions and yield of ScHerBELs expt. 5. 


The number of stages plus one is then found from 


24, 


m 
For small values of « and 4 (i.e. sharp separations) 
the total number of stages is more conveniently derived 
by an approx mation which can then be shown to be 
very accurate. (See next section.) 
It is of interest to consider the connection with the 
special solutions discussed in section 2. 


The solution m = xn corresponds with O = 0 or 


simply the t-axis. P, and P, become the intersections 


it is of course unnecessary to draw the 
From the corresponding @, values all 
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of C, and C, with the taxis. Their abscissae can be 
easily found from 


or simply 
« 
log log R, 
and in a similar way 
ty log log R, 


Hence we find for M the abscissa 4 iog R, Ry. 


which leads to 
2 log R, R, 
log 


gives ¢, 


The special solution x f, and cor- 
responds with the intersection of the C, and C, curves. 


The abscissa can be found from 


Apart from the trivial solution ¢ — 0, this has the VOL. 
solution l 

(l--a) 952 
ty blog Slog Ry Ry 
which leads to the same result 
2 log RR, 
mon log 


As an application of this method the example 
of Scnerpet already treated in the previous sec- 
tion is considered. 


R,= 0-256; R, 1-620. 


From this we get a — 0-1078 and b = 0-2038. 
The appropriate parts of the corresponding C, 

and ©, curves are drawn. The results are plotted 

in Fig. 5. As a<b we replace w by which is 


allowed in view of the symmetry of Fig. 3. 
A rather flat minimum is found. We obtain 


= 084: 0-66 and — 0-06 so that 


It is seen that the minimum number of stages plus 
one (n + m) is not substantially below that given for 
cases 3 and 4 in table 1 (viz. n — m = 14-4). 

The curves shown in Fig. 2 were derived from 
SCHEIBELs expt. |. Here the difference is somewhat 
ureater. 

4. AN APPROXIMATE SOLUTION 
-Equation (22) gives an upper estimate for the minimum 
number of stages. In order to judge its accuracy 
a lower estimate is also required. 


| 
OAT 
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For sharp separations eq.(1) and (2) may be — with 


simplified by putting m | log R, 
Ee BMI This is a lower approximation for the exact minimum 
so that The two approximations for (n m) will differ: 
» that: 
log R, . 
ky OR, | by less than 2% if the ratio of ; ye is between 
(35) << og R, 
= R, | the limits and 
, by less than 10% if this ratio is between ! and 4. 
whence, by using eq. (3) and (3a): 
So the approximation to the exact solution is 
log R log R, ertainiv hetter 
(m + n)logp (36) certainly better. 


SCHEIBELs expt. 5 (Ry — 0-256, R, — 8-28) is still 
This function is to be compared with the function within the inner limits. 


used in section 3 of which 6 special solutions were 
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NOTATION 
Symbol Definition Unit 
( 32) limension| 
R 
(eq. 32) dimensionless 
extaction factor {see part I (eq. dimensionless 
m number of stages in washing section, inel, feed stage dimensionless 
" number if stages in extraction section, incl. feed stage dimensionless 
R ratio of the amounts of solute in final extract and extraction residue [see part | eq. (15)] dimensionless 
n) [eq. (31)} 
t abscissa in Fig. 2 
ordinate in Fig. 2 [see eq. (38a) and (33b)] 
log 
og (31)] dimensionless 
log 
log 
log 
ratio of distribution coefficients dimensionless 
kK, E, 
ra 
Indices: 
1 component with lower distribution coefficient 
2 component with higher distribution coefficient 
Subscript : 
0 extraction conditions for which the total number of stages is a minimum 


= 
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Book review 


L. Polymerisationskinetik, Springer-Verlag, Hei 
delberg, 1951. 44 figures, 31 tables; VIII 287 pp.. DM 36.60 


This is the first up to date German book on the mechanism 
of polymerisation reactions. Advances in this field in the 
last half decade have been so rapid that the material is already 
seattered in the literature waiting to be collected, collated, 
and critically surveyed. This volume fulfils these aims and 
is written for physical chemists who are practising in this field 
and want to know how matters stand and also for the more 
general reader who cannot possibly keep pace with the flood 
of literature now being produced on this subject. Fortunately 
the subject can be treated in a systematic and logical manner. 
Furthermore recent developments have been to a ve ry large 
extent quantitative which makes it much more satisfactor) 


to write a monograph on a newly developing topr 


In writing a book on any branch of polymer chemistry 
it is not easy to keep strictly to the subject under review. The 
author has, however, carefully confined himself to the kinetics 
of reactions and has only made passing reference to ancillary 
subjects. There is a brief general introduction to say what a 
polymer is and to indicate what chemical kinetic analysis 
ean do for the elucidation of me« hanisms. There os too, a 
brief description of the various oNxp rimental methods for 
following polymerisation. As a necessary adjunct to these 
techniques molec ular weht mhoewsurements are nearly always 


necessary if the mechanism is to be completely understood. 


There is therefore a short section describing the usual methods. 
Space is not sufficient, however, to give practical details for 
a newcomer to the field. The real part of the book begins with 
a general chapter on addition polymerisation. The theory of 
non-stationary state kinetics is also proposed but the new 
developments in this latter field are far too recent to have 
been incorporated in the book. The major part of the book 
deals with sper ifi« examples. Here eve ry aspect of the subject 
is clearly discussed and appraised, including such developments 
as the determination of rate constants, All the developments 
in copolymerisation are thoroughly discussed and the more 
recent deve lopments on the rates of copoly mersation are not 
merely added as an addendum. The controversial and less 
quantitative subject of ionre polymerisation is also dealt with. 
In such polymerisations matters are not nearly so precise 
kinetically and correspondingly the chapter on this topic is 
brief and summarises the present state of affairs. Polyeonden 
sation kinetics of linear polymers are somewhat simpler than 
those of addition polymerisation and the now classical treat 
ment is given in detail not only for rates but also for molecular 
weight distributions. Not so much attention is given to the 
more difficult state of affairs in the synthesis of branched 
polymers. The volume is therefore an up to date survey of the 
whole field. 


ment has been omitted. 


It is remarkably complete and no new develop 
The author has welded the material 
into a very readable and satisfactory account of the subject. 


H.W. 
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